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Preparation of this document 

FAO-ICLARM Stock Assessment Tools (FiSAT) is the software that 
resulted from the merging of LFSA (Length-based Fish Stock 
Assessment) developed by FAO and the ELEFAN (Electronic Length 
Frequency Analysis) package developed at ICLARM* (International 
Center for Living Aquatic Resources Management), and the addition of 
some routines found useful for the analysis of length frequencies. 

The MS DOS version of FiSAT was converted to a Windows compatible 
version, hereafter referred to as FiSAT II. This work was carried out 
between 2000 and 2002 through the European Union support project, 
Fisheries Information and Analysis System (FIAS). 

Development of the FISAT II software was carried out by Felimon 
Gayanilo. The new software was thoroughly tested during an FAO 
working group meeting at Cidade da Praia, Cape Verde, in August 2003, 
under the leadership of Ana Maria Caramelo and Merete Tandstad, from 
the FAO Fisheries Department, producing a detailed report of the major 
errors then detected in the software. 

The errors and bugs in the software detected during the workshop, as 
well as those that were discovered during the later testing phase were 
corrected by Pedro de Barros, Professor at the University of Algarve. 

This User’s Guide is an update of the previous version (Gayanilo, Sparre 
and Pauly, 1996). In addition to the printed document, a compact disk is 
provided with the software. 
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ABSTRACT 

FiSAT II is a program package consisting of methodologies for 
use with computers, enabling users to formulate some 
management options for fisheries, especially in data-sparse, 
tropical contexts. 

The FiSAT II was developed for computers running on Microsoft 
Windows operating systems. The new version utilizes the 
standard Windows graphic user interface. 

FiSAT II was developed mainly for the analysis of length- 
frequency data, but also enables related analyses, of size-at- 
age, catch-at-age, selection and others. 
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Chapter 1. Getting started 


What you will learn from this chapter 

In this chapter, you will be introduced to terminologies 
used in this guide, shown how to install the package, and 
understand the objects used in the graphic user interfaces of 
FiSAT II. At the end of the chapter, you will be introduced 
to the revised menus of FiSAT. 


About the package 

FiSAT (FAO-ICLARM Stock Assessment Tools) is the 
product that resulted from merging the LFSA (Length- 
based Fish Stock Assessment) developed at FAO (Sparre, 
1987) with the Compleat ELEFAN (Electronic Length 
Frequency Analysis) package developed at ICLARM* 
(International Center for Living Aquatic Resources 
Management) (Gayanilo, Soriano and Pauly, 1989). All 
routines in the earlier packages are included in FiSAT II. 

Descriptions of many models used in this package may also 
be found in Pauly (1984a), Sparre and Venema (1992) and 
the FiSAT reference manual (Gayanilo and Pauly, 1998). 


System requirements 

FiSAT II, which is distributed on a compact disk, can be 
installed and made executable on many computer platforms 
or computing environments that run on the Microsoft 
Windows NT/2000/98/95 operating system. The minimum 
system requirements are as those identified when installing 
Microsoft Windows operating systems: 

* WorldFish Center. 
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• Microsoft Windows 95/98/2000, 

• 64MB of RAM, 

• 1 024 x 768 high-resolution monitor, and 

• at least 5 MB free disk space in the Windows 
directory and another 6 megabytes of free disk 
space for the destination address. 


Installing and starting FiSAT II 

The CD-ROM that comes with the package contains the 
setup.exe. Run the SETUP.EXE found in the root directory 
of the CD-ROM. The installation wizard will require user 
responses to prompts (Fig. 1.1). 


DO NOT use other methods for installing the package to a 
hard disk. The SETUP program does not only transfer the 
files from the source to the destination disk but also 
initializes other files to be used by FiSAT II. 


The installation routine will also install sample data files 
that can be used and an item (icon and title) in the 
Windows Programs menu. Selecting the item activates 
FiSAT II. 
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Welcome to the FiSAI II Setup 
Wizard 

This will install FiSAT II (Vet 1 .1 .0) on your computer. 

It is recommended that you close al othei applications before 
continuing 

Cick Next to continue, or Cancel to exit Setup. 


Next > I Cancel 


Fig. 1.1. A graphic user interface of SETUP.EXE when installing the FiSAT II package, 
the directory should be: \Program Files\FiSAT II. 


It is important to note that the first display screen (hereafter 
referred to “Splash Screen”; Fig. 1 .2) contains the version 
installed. This number (also available in the About... 
form), is used in reporting programming defects or in 
accessing technical support from the authors. 



Fig. 1.2. FiSAT II splash screen showing the version number of the software 
installed. 
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The graphic user interface 

FiSAT II utilizes standard icons and objects found in most 
commercial software packages that run in the Microsoft 
Windows environment to make it as user-friendly as 
possible (please refer to Microsoft Windows documentation 
for details). 

The main user interface 

The main features of the FiSAT II main user interface (Fig. 
1.3) that may differ from other commercially available 
software packages are: 


1. Windows Caption: The content of the caption also 
displays the version number and in addition, the caption 
of other windows loaded into memory if windows is 
maximized to occupy the whole working area. 

2. Menu Bar: This bar contains the various menus 
available to the user in accessing commands (see 
further below). 

3. Tool Bar: This bar contains command buttons to 
commonly used functions. Double-clicking the Tool 
Bar will open a Window that will allow the user to 
customize the contents of the Tool Bar. 

4. Status Bar: This bar contains immediate help messages 
in the left panel (see also Help Messages below), the 
system date and time in the right panels used in printing 
of reports and time stamping of files. 
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Fig. 1.3. The main window of FiSAT II. Please refer to Microsoft Windows 
documentation for detailed explanation of other standard objects found in a window. 


User inputs to FiSAT II 

To facilitate the encoding of data, input boxes (Fig. 1 .4) are 
color-coded. Input boxes with 

• Yellow background are required fields, 

• White background indicates fields that are optional, 
although it is advisable to encode data whenever 
available to prevent data losses, and 

• Gray background indicates fields that cannot be 
edited by the user. 

FiSAT II uses standard fonts for the user interfaces. 
However, FiSAT II automatically converts font colors to 
red if the value was computed by the software. Also, 
grayed labels on command buttons of filed names indicates 
that the object is presently not available to the user. These 
are enabled only if the pre-requisites for the function to be 
available are met satisfactorily. 
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Yellow background User 
required to input valid 
values 

Gray background Values 
cannot be edited 

White background: Input 

to these fields are optional 

Red Font Value computed 
by FiSAT II 

Grayed Labels: Function is 
disabled 


Bint- 


0o« 


Fig. 1 .4. Typical user interface showing fields with yellow, white and gray backgrounds. Also, 
note the use of the red font to indicate fields computed by FiSAT II. Also, note that in some 
cases, command are presently not available to the user. Selecting points directly from a graph 
are available. 


There are a number of routines in FiSAT II requiring the 
user to select a point in the graph. Move the mouse pointer 
to a point to select. The standard cursor icon will change to 
a cross (Fig. 1.4) if a point can be selected. Clicking the 
mouse will select the point. If a point is selected, the yellow 
circle will be shaded with black. 


Using the menu bar 

The Menu Bar of FiSAT II contains the six other menus. 
Clicking on the menu will dropdown options. For example, 
diking on the “File” option will drop the commands 
available to the user (Fig. 1.5). Note that some of the 
selections are disabled (colored grayed) indicating that the 
function(s) are currently not available. If submenus are 
available, an arrow head is displayed in the right comer of 
the menu item. Clicking the items with submenus will 
dropdown the contents of the submenu. 
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•jFiSAT II (Version I. I.O)... 


File Ed* Assess Support Window Hdp 


^iSJx] 



Fig. 1.5. Menu system of FiSAT II. Note that some functions are disabled if it does not 
apply for the given situation (in this case, no file is loaded). Also, note the presence of 
the arrowhead (in red circle) if a submenu is available. 


Alternatively, pressing the <ALT> key together with the 
underlined character of the menu item access the command. 
For example, pressing <ALT> + <P> will open the Print 
dialog box. Another option is using the shortcut keys (not 
available to all commands). In Fig. 1.5, pressing 
<CTRL>+<P> will also open the Print dialog box. 
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Chapter 2. HELP facilities 


What you will learn from this chapter 

In this short chapter, you will be introduced to the types of 
help messages incorporated in FiSAT II, to assist you in 
using its various routines. However, these help messages 
are not intended to replace this User's Guide. 


Types of HELP messages 

FiSAT II has three types of help messages, classified as 
either immediate, warning or descriptive messages. 

• Immediate help messages : These automatically appear 
in the Status Bar (i.e. at the bottom of the screen), when 
certain actions are to be performed. They consist of one 
line of information (less than 66 characters), generally to 
remind users of the limits for the current entry field or what 
options are available. 

• Warning messages : These are messages that appear at 
the bottom of the screen to warn users about the value 
entered in the current field. In most cases, FiSAT II may 
require the user to provide another value before a function 
is performed. 

• Descriptive help messages : These describe the routines 
under each item of a given menu, or provide a list of keys 
in a given routine that can be used but which are not shown 
on the display screen. 
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Accessing HELP 

The Help Menu of FISAT II provides three ways of getting 
help as the need arises. 

On-line help messages 

Accessing the “Contents” menu item in the Help menu or 
pressing <F1> will open the on-line user guide (Fig. 2.1). 
The contents of this manual were electronically converted 
to facilitate the access of descriptive help messages. 

<9 u- 



Fig. 2.1. FiSAT It on-line help messages. This is the 
electronic version of this document. 


This standard Microsoft Windows type interface provides 
facilities to view the pages through the Contents list (as in 
Fig. 2.1) or specific topics can be located via an index field 
or free text searches. Please refer to Microsoft Windows 
references for details on how to use this interface. 
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FiSAT II on the web 

A FiSAT II website is maintained by FAO to provide a 
venue for distributing updates to the software and other 
announcements related to the application of the software 
(Fig. 2.2). 



Fig. 2.2. FiSAT II home page. The content of the pages varies as the need arises. 
The current URL is http://www.fao.org/fi/statist/fisoftyfisat/index.html. 


The content of the pages changes as the need arises. This 
website will contain pages to (i) allow users to quickly 
assess the new package, (ii) online user guide with the 
same content as the online help messages accessed through 
the Help Menu, (iii) user registration to allow FAO to 
update users on new developments, (iv) pages to download 
new updates as they become available, and pages to answer 
frequently asked questions (FAQ). 

Users are encouraged to register their names and to 
regularly visit the website to keep phase with the 
developments of the software. 
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Technical support 

In cases where the online help messages and the FiSAT II 
website are not sufficient to answer a specific query, or in 
reporting programming defects, FiSAT II users may 
forward their queries to FAO. To facilitate the services, 
always include the following: 

1 . The FiSAT II version installed; 

2. The Microsoft Windows operating system; and 

3. Step-by-step procedures to repeat the problem 
(if applicable). 

In rare cases, FAO may request users to access and e-mail 
system information. These files may be accessed from the 
“About FiSAT II” form accessible through the Help Menu 
(Fig. 2.3). 

The FAO technical staff will provide a step-by-step 
procedure on what needs to be retrieved and mailed to 
resolve a technical problem. Please refer to Microsoft 
Windows reference manuals for details on the use of the 
interface that displays system information. 



*J 


FiSAT II 

FAO-iCLARM Fish Stock Assessment Tools 


Warning All rights reserved No pari of this product 
map be reproduced or transmitted m any form without 
prior permission ot the copyright owner Application tor 
such permission, with a statement of the purpose and 
extent ot the reproduction, should be addiessed to the 
Duectoi. Publications Division. Food and Agncdture 
Organization ot the United Nations, Vialle dele T etme 
di Caracalta, 00100 Rome, Italy 
Copyright *2000-2005. FAO. Rome. Italy (Vetsionl l O) 


[ "OK I 

System Info... 


Fig. 2.3. The "About FiSAT II" form. The system information may also be 
accessed from this form. 
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Chapter 3. FILE menu 


What you will learn from this chapter 

In this chapter, you will be introduced to the different types 
of data that can be created and used in FiSAT II, how these 
data are encoded, edited, imported, exported and options 
available for manipulating data files. 

Types of files 

FiSAT II supports several types of data differentiated 
largely by their filename. Table 3.1. is a summary of the 
different file-types that can be created in FiSAT II. 


Table 3.1 . Summary of file-types supported in FiSAT II and their corresponding 
file extensions (see details of each file type below) 


File -Type 

File Extension 

Grouped length frequencies 

Mfq 

Ungrouped length frequencies 


Variable grouping of length frequencies 


Length-at-age 

Maa 

Growth increments 

‘•gin 

Catch-at-age 

*.caa 

Mean and s.d. 

*.msd 

Probabilities 

*.prb 

Value by size 

*.val 

XV regression file 

*.xyf 

Species table 

*,spc 

Fleet table 

*.flt 

Multispecies/multifleet scenario 

‘.see 
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Note that, although FiSAT supports a number of data types, 
most of the routines available in FiSAT either require 
length-frequency data or data derived from length- 
frequencies. The different data files that can be created and 
stored using FiSAT are: 

Time series of length frequencies 

Data files of this form can either be grouped by constant 
class size (i.e. class interval), grouped by variable class 
sizes (common for data from commercial sampling of 
shrimp) or "ungrouped frequencies", i.e. individual 
measurements requiring subsequent grouping. 

Length frequencies grouped by constant class size 

Length frequencies of this type are characterized by lower 
and upper class limits, and a "mid-length", the mean of the 
limits (Fig. 3.1). FiSAT II does not require inputs for the 
upper and lower class limits but rather uses the class size 
(or intervals) and the smallest mid-length, which enables 
the mid-lengths of the subsequent classes to be computed. 

Unlike other data types, this data requires inputs to sample 
sub-headers (Table 3.2) that include, amongst other things, 
the sampling date in the format (dd/mm/yyyy). The data 
file will not be properly saved if the sampling date is not 
encoded properly. Samples with invalid or missing 
sampling dates will be automatically deleted. A command 
button is available in the user interface (Fig. 3.1) to open 
the sample sub-header. Alternatively, double-clicking on 
the column or clicking on the column header opens the 
sample sub-header. 


Copyrighted material 



I 


14 


Lr:..Pro(|r«m F«n\riSAT n ;Dal«.,SnmpfcJfq 


^131*1 


Speaks name |ljpenrm moitccenw 

R<nrarit 


Imported fionri F (SAT 
(DOS vei J: Thu 0 0 
lest dota 

Oihet (dsruhej |RAGAY GULF 1980 

linfolmMtucerort | OT zi Type ol meuuma* |Tc«llerWhlTL) _l| 

Srwolleil mdfcnglb [5 5 □on rterval [t 


r 


, e 

Edt Sarple Header | 

AddSampfc 

Delete “ample 

I P Duplay decr-i^ ponton 



ML 

20/1/1985 

2/3/1980 

21/5/1980 23/7/1980 

27/9/1980 20/11/1980 - 


5.5 

'I 








6.5 

1 





3 


75 

4 

1 




9 


85 

22 


3 



9 

95 

213 



6 

19 

4 

10.5 

105 

77 

26 

26 

13 

108 


11.5 

103 

198 

423 

61 

37 

245 


12.5 

8b 

62 

415 

148 

44 

235 


13.5 

68 

58 

289 

138 

115 

133 


145 

53 

24 

185 

75 

80 

107 


155 

15 


40 

32 

38 

30 

165 

3 


36 

21 

8 

10 

17.5 



3 



4 

18.5 















< 

u 






lC 


Fig. 3.1. User interface for recording grouped length frequencies. Note the command 
button to edit the sample sub-header. Alternatively, double-clicking on the column opens 
the form to encode the parameters for the sub-header. 


The required inputs are: 

• Date of sampling (DD/MM/YYYY) 

• Smallest class mid-length (from 0.25 to 500) 

• Class size (from 0.25 to 500) 

• Type of measurement 

• Frequencies (0 to 999,999) for each class. 


WARNING! Although it is possible to copy contents from other 
Microsoft Windows software through the Microsoft clipboard 
(i.e. Copy->Paste procedure), such as Microsoft Excel, care 
should be taken to check for the input to sampling dates. 
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Table 3.2. Sub-header for size (length or weight) frequency data grouped by constant class 
interval. 


Field 

Type 

Remarks 

Sample No. 

N 

The values for this field are automatically 
provided by FiSAT and cannot be edited. 

Day of sampling 

N 

The upper limit is dependent on the month 
and year. A default DAY (equal to 15) is 
provided for convenience. 

Month of sampling 

N 

The default value is the previous month + 1 . 

Year of sampling 

N 

Year of sampling, from 1900 to 2100; the 
default for this field is the year of the 
previous sample. 

Sex* 

C 

Valid entries to this field are: female, male 
or undefined/both. Only the first letter of 
the choice needs to be entered; FiSAT will 
complete the entry when <ENTER> is 
pressed. 

Depth (m) 

C 

This field should indicate the (mean) depth 
at which the sample was obtained. 

Area 

c 

This indicates the area from which a sample 
originates. 

Latitude 

c 

The input to this field, indicating the 
location of the sample will not be validated. 
However, it is recommended that a standard 
format be used, e.g., XDD:MM:SS, where 
X refers cither to N (north) or S (south), DD 
refers to degrees oflatitudc, MM to minutes 
and SS to seconds (example: S35:30:00 ) 

Longitude 

c 

As above. Thus, YDDD:MM:SS, 
where Y is cither E (east) or W (west), 
DDD refers to degrees of longitude, MM to 
minutes and SS to seconds (example: 

W12: 10:30) 

Gear type 

c 

This field identifies the gear type used for 
sampling (example: Trawl). 

Gear characteristics 

c 

This field should describe specific 
characteristics of the gear used 
(example: 2.5 cm mesh) 

Vessel 

c 

This field should contain the vessel name or 
type. 


+ C - character field; N - numeric field; Ch - limited choice field 
♦The fields from Sex to Remarks arc optional fields 
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Table 3.2 (Cont). Sub-header for size (length or weight) frequency data grouped by 
constant class interval. 


Field 

Type* 

Remarks 

Time of day 

C 

This field refers, on a 24h basis, to the time 
the collection started or ended in hours and 
minutes, (example: 14:30). 

Bottom type 

C 

This field indicates the bottom structure at 
the sampling site (example: sandy). 

Type of catch 

C 

This field may be used to describe the data 
at hand, e.g., commercial or research data. 

Effort 

N 

This field should contain a measure of the 
effort expended to obtain the sample in 
question (e.g., no. of hauls, no. of hours, no. 
of boats, etc.). This field will be used when 
the data arc to be re-expressed as catch per 
unit of effort (C/0; thus it is important that 
the unit used is standardized. 

Wt. of the sample 

N 

This field should indicate the weight of the 
sample, as taken out of the total catch. This 
implies that the unit used in measuring 
weight should be consistent for all samples. 

Total wt. of species 

N 

This field should indicate the weight of the 
species as a subset of the sample collected. 
It is important to note that the value 
provided should be less than the WT. OF 
THE SAMPLE and that the unit used in 
measuring weights should be the same 
throughout. 

Total wt. of catch 

N 

This field is for a value representing the 
total multi-species catch. As in the two 
previous fields, the unit used in measuring 
weight should be the same. 

a 

N 

The coefficient "a" of the length-weight 
relationship (W=aL* 5 ). The default value is 
the one used in the previous sample. 

b 

N 

The coefficient "b" of the length-weight 
relationship complementing the previous 
field. The default value is 
the one used in the previous sample. 

Remarks 

C 

Other information not included in any of 
the above fields. 
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Ungrouped length frequencies 

This form of length frequencies is characterized by only 
one column of inputs (Fig. 3.2), the length of individual 
specimens. As with the previous data form, data of this 
form cannot be used for analysis in FiSAT II unless 
regrouped with constant class size. This is provided as an 
option in FiSAT II when saving this type of file. 


C:\PriMjrant riles \Fi5AT 1 1 \U . 1 1 a \Samp le.ulf 




-ini xi 


Headers 


Species name 
Other identifier 
Sarrplins date 


lllpeneus moluccensis 
IRaflay Guf 1980. Phiippines 


15 Julho 1981 

Unit of measurement Type of measurement |lotal length |TL) 




rrequana 


Ob* 

Length (cm| 

1 

iTool _ 

2 

sou 

3 

9.00 

4 

12 00 

5 

12.00 

G 

12.00 

7 

is on 

8 

7 00 

9 

9 00 

10 

IG.OO 

11 

14.00 

12 

14 on 

<1 




Remarks 


1 his is only a sample data. 


Fig. 3.2. User interface in encoding ungrouped length frequencies. This data type may 
also be saved as grouped length frequencies. 


Copyrighted material 




18 


The required inputs are: 

• Date of sampling (DD/MM/YYYY); 

• Measured lengths (0.25 to 999) 

When saving to a form with a constant class size, the 
required inputs are: 

• Class interval to use (0.25 to 500); 

• Smallest class mid-length to use (0.25 to 500); 

• Filename to store the new file (if the filename 
already exists, the data will be appended to the 

existing file). 


Data grouped bv variable class sizes 

This form of length frequency data is characterized by 
explicit lower and upper limits for each class (Fig. 3.3). 
The mid-length cannot be used in this case. Files in this 
form can be stored, but they cannot be used within FiSAT 
II unless regrouped with constant class size. 
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* l >D\ Program r ilcs r iS AT II • Doto S ample. vll 


- '□! x| 


Species nor® 

Other identifier 

Sampling date 
Unit of measurement 


|Upeneus molucxensis 
|Ragay GuB. Philppines 1 960 
1 15 JJho 1980 _* *| 

foe Type of measurement | Total length |TL] 


3 


I - Dcplay decimal position 



Remarks 

This is a test data few variable 
length groups. 


Fig. 3.3. User interface to encode length frequencies with variable length classes. 


The required inputs are: 

• Date of sampling ( DD/MM/YYYY); 

• All lower and upper class limits (0.25 to 999); 

• Frequencies (0 to 999,999) for each class. 

When saving to a form with a constant class size, the 
required inputs are: 

• Class size (0.25 to 500); 

• Smallest mid-length to use (0.25 to 500); 

• Filename to store the new file (if the filename 
already exists, the data will be appended to the 

existing file). 
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Growth increments 

Growth increments are usually referred to as "mark- 
recapture" data (because mark-recapture studies are a major 
source of growth increment data). Pairs of mean lengths 
that are linked in modal progression analysis may be saved 
in FiSAT II and analysed as growth increments. 

This data type is characterized by four columns of entries 
(Fig. 3.4); the initial date of observation (relative to 
marking), the initial size, the date of recapture or date of 
second reading of the same specimen and the size at 
recapture or second reading. FiSAT II will accept lengths at 
recapture implying zero or negative growth, but will reject 
a date pertaining to a second reading that is earlier than the 
initial date or date of first reading. 


' C:\Progr am nlef\FISAT II'. Data' Sample t.gln 


^JOjxJ 


Species name 
Olhei iderliliw 


(Large lish 


jNewmajoi exercise (bavffei) 

Unit of length measurement »] Type of measmement |Olhers/Unknown 3 


Frequencies 

(NOTE Yellow cdumn refers lo tagging dale and unshaded le whee. tefere todala 
5* receive I 


Obs 

Dur 

Month 

Ye* 

Length 

Ddv» 

Month 

Yeai 

1 *L 

1 

is| 

i 

1965 

14.98 

15 

2 

1985 

— 

2 

15 

2 

1985 

IB 08 

15 

3 

1985 


3 

IS 

3 

1985 

17 08 

15 

4 

1985 


4 

15 

4 

1985 

18.34 

15 

5 

1985 


5 

15 

5 

1985 

18 78 

15 

6 

1985 


6 

15 

e 

1985 

70 13 

15 

7 

1905 


7 

15 

i 

1985 

0.84 

15 

2 

1985 


a 

15 

2 

1985 

10.29 

15 

3 

1985 


9 

15 

3 

1985 

11 37 

15 

4 

1985 


hi 

15 

4 

1985 

1232 

15 

5 

1985 


it 

15 

5 

1985 

12.92 

1b 

8 

1985 


12 

15 

6 

1985 

14 58 

15 

7 

1985 


iL 







► 

E 


Remarks 


Impeded bom 
FiSAT (DOS 
ver J, Data 
computed from 
length 

frequencies: 
Generated 
from hiking of 
me are 


Fig. 3.4. User interface in encoding growth increments. The columns with a dark yellow 
background refer to tagging data and those with a white background are recaptured data. 
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The required inputs are: 

• Initial date or date of first reading 
(DD/MM/YYYY); 

• Corresponding length (0.25 to 999); 

• Final date or date of second reading (DD/MM/YYYY ); 

• Corresponding length (0.25 to 999). 


Length-at-age 

The age of a fish may be obtained by directly reading 
otoliths, scales or other hard parts. Pairs of observations, 
(age, length) collected for individual fish, form length-at- 
age data. Series of mean lengths, as identified by modal 
class progression analysis, can also be treated as length-at- 
age data, although they actually represent lengths at 
relative age (precluding the direct estimation of t 0 ). 

This type of data is characterized by two column entries; 
the ages (or relative ages if derived from modal progression 
analysis) and the corresponding length measurements 
(Fig. 3.5). The ages can either be in days, months or years 
and may be any values from 0 to 999. The length values 
can only be from 0.25 to 999, and the units must either be 
in cm, mm or inches. FiSAT II provides a third column 
entry for a weighting factor (WF) for each observation. 
These weighting factors (default of 1.0; acceptable range 
0.01-100) can be used to emphasize or de-emphasize 
observations during data analysis. 
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> C:\Progr f Hm\FISAI n.Data\Sample.lao 


-ini x[ 


Species name 
Other idertfe 
Sampling date 


|Laige species 
|Sample data ffcs 


[To Julho 2000 ^ | 

Unit ot length measurement [cm ~3 T ype ol measurement | Total length HL) ^ 
Unit ot age measurement [years » | 


Qbt 

Age 

Length 

WF. ± 

1 

0.1081 

2.000 

1.00 _ 

2 

0.108 

2 000 

1.00 

3 

0 185 

4 000 

1 00 

4 

0 270 

5 360 

1 00 

5 

0 351 

6 440 

1 00 

6 

0.437 

7 530 

1 00 

7 

0 518 

8 930 

1 00 

8 

0 603 

9 930 

1 00 

9 

0 689 

11 100 

1 00 

10 

0.771 

12 050 

1 00 

11 

0.939 

13 860 

1 00 

12 

1 024 

14 860 

1 00 

13 

1.108 

IB 020 

1 00 

14 

1.184 

17 170 

1 00 

<1 



ill 


Remarks 



Fig. 3.5. User interface in encoding length-at-age data in FiSAT II. 


The required inputs are: 

• Ages (0.01 to 999); 

• Lengths (0.25 to 999); 

• Weighting factors (default is 1.0; 0.01 to 100) 


Catch in numbers by age 


This data type (also referred to as catch-at-age. Fig. 3.6) 
can be used only for one routine, i.e. age-structured VPA. 
The required inputs are similar to length-at-age data. The 
required input consists of two columns: (i) the ages 
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(starting with the youngest group) and (ii) the catches in 
numbers. 


•^C:\P»ogr«m nir**n«5AT IT\fca*«v 5nrnplc.r «n 




[Sartirw pfchardu* 


| Caught off Morocco 


Speews name 
Other denser 
Ur\< ot c«*ch rrreasutecneri |njmbeis 
Unit ot age n*a«j'emew I years 


■a 


0b» 

Aye 

Latch 

WF. * 

1 

i «»i 

17538 000 

1 00_ 

2 

1 250 

7859 000 

1 00 

3 

i sun 

5225 000 

1 00 

4 

1 750 

34276 000 

1 00 

5 

2 0(1(1 

118049 000 

1 00 

6 

2 250 

39819 000 

i no 

7 

2500 

37976 000 

1 00 

0 

2 750 

434354 000 

1 00 

9 

3 000 

735246 000 

1 00 

10 

3 250 

383842 000 

1.00 

ii 

3 500 

106583 000 

1.00 

12 

3 750 

319512 000 

1.00 

13 

4 ono 

514256 000 

1 00 

14 

4 250 

33191 000 

1.00 

15 

4 500 

66207 (100 

1.00 

16 

4 750 

139836 000 

1 00 

17 

5 000 

15624 000 

1.00 * 

* 1 

>r 


Remark* 

Sample data from Pauly 1 964 

IM-08: Fc-0.81 


Fig. 3.6. User interface to encode catch-at-age data in FiSAT II. 


The required inputs are: 

• Age (0.01 to 999); 

• Catches, in numbers ( 1 to 999,999). 


Series of mean lengths and standard deviations 

In most cases, the data in this file type will be products of 
the decomposition of normal distributions using, e.g.. 
Bhattacharya's method or Hasselblad's NORMSEP as 
implemented in FiSAT II, which automatically generate 
such files (see next chapter). However, creation of such 
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files from external sources for subsequent analysis by 
FiSAT II is possible. 

The file is characterized by three columns of entries; (i) the 
date of sampling, (ii) the mean length, and (iii) the standard 
deviation (Fig. 3.7). 


f : Program nlr* ,rt5AT IT, Data'. Sampir 


jJBJ * 


Sp*CH»S rMTTM 
Ofr*f identfier 


|l m j* Inh 


|Ne*v nr * a It>«*lar| 

Un* d torglh rnMtuiMncrrf 


Ob* 

Dot 

Mon* 

Year 

Moan 

• d 

• 

ZZED 

1 

1905 

0 04 

1 09000 _ 

2 

15 

1 

1905 

14.90 

1 82000 

3 

15 

1 

1‘J05 

20 95 

7 ta>TMJU 

4 

15 

2 

1985 

2.00 

0 58000 

5 

15 

2 

1905 

10.29 

1 3MI0U 

6 

15 

? 

1905 

10 00 

1 04000 

7 

15 

3 

1985 

4 (111 

0 48000 

8 

15 

3 

1985 

11 37 

1 54000 

9 

15 

3 

19U5 

17.08 

1 90000 

10 

15 

4 

1985 

537 

0 88000 

11 

15 

4 

1985 

12.32 

1 65000 

1? 

15 

4 

1985 

18 34 

1 81000 

13 

15 

5 

1905 

G4S 

0 70000 

14 

15 

5 

1985 

12 92 

1 31000 

15 

15 

5 

1985 

18 /0 

7 290110 

IS 

15 

6 

1985 

7.54 

0.93000 

17 

15 

6 

1985 

14 58 

1 99000 

le 

1 

15 

6 

1985 

20 13 

2 29000 „ 

it 


Imported tram FiS AT |DOS 
I Data computed bom 

length Itequenoe* 


Fig. 3.7 User interface in encoding means and standard deviation file. 


The required inputs are: 

• Date of observation (DD/MM/YYYY); 

• Mean length (0.25 to 999); 

• Standard deviation. 


Probabilities of capture 

One common cause of bias in length-frequency data is the 
selectivity of the gear(s) used to obtain the samples. This 
bias can be in part overcome by correcting size frequency 
samples with the probabilities of capture stored in a 
"probabilities of capture file", and generated either through 
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a selection experiment, or a length-converted catch curve 
(see next chapter). 


This type of data file (Fig. 3.8) is characterized by two 
columns of entries; (i) the class mid-length (0.25 to 999) 
and (ii) the probability of capture which may range from 0 
to a maximum of 1. Entries less than 0.001 will be 
registered as 0. Note that a file in which all values are equal 
to 0 or 1 will not be stored. 


USES 


,«smss!iafl3 


^injxj 


Spaces name 
Other idenrjfiei 


|fjpeneu& moluccercts 


|RAGAY GULF I960 
Unit of length r»e«iremen( "on »] 


pFia 


nb< 

1 

n- ■ . .i 

a n 

PiaheUlilr 

0 2405 7 
0 35676 

iTbul 

2 

6.50 

3 

7.50 

0.48035 

4 

B.50 

0.55026 

5 

9.50 

0.73198 

6 

10 50 

0 78741 

7 

11 50 

0 91668 

a 

12.50 

1.UU0U0 


id 




Estirrvsied from Lengtv 
converted Catch Curve 


Fig. 3.8. User interface in encoding probabilities in FiSAT II. 


The required inputs are; 

• Smallest class mid-length (0.25 to 500); 

• Class size (0.25 to 500); 

• Probability of capture (default is 1 .0; 0 to 1.0) 
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Table of values 

The commercial (e.g., market price) or biological value 
(e.g., in terms of fecundity) of fish often varies with size. 
This can be accommodated, for various assessments by 
using size-specific values as multipliers. In FISAT II, the 
values are used as a multiplier in the Thompson and Bell 
Yield-Stock Prediction routine. 


This data type (Fig. 3.9) is characterized by two columns of 
data: (i) Class midlength, and (ii) value for each group. 


C:\Progrwi Files \F BAT II Data 


•.la*! 


Specie* name 
Olhei idenlifiei 


[lipeneu* rrarfucerai* 
(RAGAY QUIP 


Unit of 'neasurefnerr [cm 3 Cutency |0BD 


Obi 

Mkflonylh 

Value 


1 

550 

1 00 



2 

C 50 

1 00 



3 

750 

1.00 



4 

8 50 

1 (10 



5 

9 50 

1 00 



6 

10 50 

1.00 



7 

11.50 

1.00 



8 

1260 

1 00 



9 

1350 

1 00 



10 

1450 

1.00 



11 

15 50 

7 00 



12 

1650 

200 



13 

17 50 

2.00 



14 

« | 

19.50 r - 

2.UU 

_£j 

t 




Rp'JprctitPfP 


Fig. 3.9. User interface for encoding table of values. 


The required inputs are: 

• Date of sampling (DD/MM/YYYY); 

• Lower and upper class limits (0.25 to 999); 

• Value (0 to 999) for each class. 
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Two-column regression files 

Linear regression analysis is an important statistical tool for 
fish stock assessment. One common application is to derive 
the relationship between the (log- transformed) lengths and 
weights of fish. 


This data type (Fig. 3.10) is characterized by two columns 
of entries, the X-values and the Y-values, both of which 
may range from 0 to 999. In FiSAT II, only the regression 
analysis routine can use this type of data. It is also the only 
routine whose headers are used only for labelling graphs. 


■ -C:\Proqriw Fllet\FISAI MDat.i Sonnilr.Kyf 


-ini xi 


Graph T<le l^geLength Reiatorei^T 
X-am Labai |Age(ye»*i 
Y-aw Label | Length (cm) 


Ob* 

X 

Y 

± 

1 

oImI 

2 000 


2 

0185 

4 000 


3 

0 270 

5 360 


4 

0 351 

6 440 


5 

0 437 

7 530 


6 

0518 

8 930 


7 

0 603 

9 930 


8 

0 689 

11 100 


9 

0 771 

12 050 


10 

0 939 

13 860 


11 

1 024 

14 860 


12 

1 108 

16 020 


13 

1.184 

17.170 


14 

1 269 

17 850 

* 



* 



Fig. 3.10. User interface for encoding two-column data for regression analysis. 


The required inputs are: 

• X-axis and Y-axis labels; 

• X-values (0 to 999); 

• Y-values (0 to 999). 
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Species table 


The Thompson and Bell Yield-Stock Prediction routine 
requires population parameters and importantly, an array of 
fishing mortality (F). In FiSAT II, these data are results 
when executing length-structured VPA. To simplify the 
procedures in creating a scenario for the Thompson and 
Bell Yield-Stock Prediction routine, these results are saved 
as a Species Table (Fig. 3.11 and Fig. 3.12). Note, 
however, that these data may also be manually encoded 
and/or edited by the user. 


pm.*: ; j: n 

| £naj«] 

Species Name 
Remarks 


[i pennus mtijccenas 


F lie generated Itom len£*> strxiued VPA. 


Fie Delab 


SmaHesi Cbts Lower Lrrw 
Clast toe 

bUrt** n! fertfh 91014* 


PcfxJaton P jwr«l<r . 


I eanrd Fishing Martaby (Ft) | 
Nato-dMotaby |M| 

Aiyrrtr c*etc length [Loo] 
VBGFGiovdi Constant (K) f 
Reciutmenl lolheFijheiy | 


LengttvWo^l PeWiorsbp 

'a* Coallocrf f~ 


b'CoalficKnt 


*J 


!] 


Fig. 3.1 1. User interface for encoding population parameters when 
creating/editing Species Tables. 
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•C:\Progiom rile5\riSAT II\Dota\5amplc5PC 


Head® 

Anays] 




Obs 

From 

To 

F-array 

Value 

1 

5^0| 

6.50 

0 0006 

1.00 

2 

6 50 

750 

00026 

1 DO 

3 

7.50 

8.50 

0 0126 

too 

4 

850 

9.50 

0.0282 

too 

5 

9.50 

10.50 

0.2150 

too 

6 

1050 

11.50 

03751 

too 

7 

11 50 

1250 

1 3314 

1 no 

8 

1250 

13.50 

1 8474 

too 

9 

13.50 

1450 

2.4383 

too 

10 

14.50 

1550 

30168 

too 

11 

15.50 

16.50 

21531 

1UU 

12 

16.50 

17.50 

2.9260 

too 

13 

17 50 

18.50 

0 0000 

ooo 

14 

18 50 

1950 

00000 

0 00 


Fig. 3.12. Portion of the user interface (second tab as shown 
in Fig. 3.1 1) to record F-array and table of values. 


The required inputs are: 

• Population parameters (class details, F,, M, L», K, 
number of recruits, and coefficients of the length- 
weight relationship [8a, b]); 

• F-array; 

• Table of values. 


Fleet table 

In simulating a multi-fleet scenario when executing the 
Thompson and Bell Yield Stock routine, a fleet indicator 
file is required to differentiate the catches recorded by 
various fleets. In FiSAT 11, a Fleet Table can be created for 
the purpose (Fig. 3.13). 
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f: .Ptrwiram HrV.riSAT II n»t/\ SwnpliM* 


Hsn« d (He FVeH ba*»th :>*wl fimi 




The fleet i* cofrpcied of m fleet* of cafctieii and o toader. The ttarded o*e*h w «d n 
Gem and no additional rratoMb « gadget are ntlalied 


d 


Fig. 3.13. User interface to encode the descriptive parameters defining a fleet. 


The required inputs are: 

• Fleet name; 

♦ Fleet description. 


Multispecies/multifleet scenario 

The file that sets the relationship between the exploited 
species and fishing fleet is the scenario file (Fig. 3.14 and 
Fig. 3.15). To compute the contribution of each fleet to the 
species’ F-array, catches by fleet are required. To allow the 
simulation on the change of mesh sizes, selection 
parameters of each fleet to a specific species also have to 
be encoded. 
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*J 


&»'«d 1 £dechon Oata and Catch indcet | 



Save Fife { OpcnFIc CieafcnHcfr | 


Fig. 3.14. Header records of the Multi-species/multi-fleet scenario file identifying 
the list of species and fleets to be included in the simulation. 
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2 
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3 

7 so 

ooo 

4 
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ooo 

5 
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ooo 

6 

IOlSO 

ooc 

7 
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ooc 

6 

1250 

000 

9 

1150 

ooo 

10 

14.50 

ooo 

11 

1550 

ooo 
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Whan me mean annual taflehei ‘a a 
apeeiat sy a fee* ye not e,-«aa&e the 
eaten *efce» >e ret«y* 
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e*r>vro a toe raw at eMNNdto o*»ar 
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each*! 
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Fig. 3.15. Tab to encode the catches and selection parameters for a species by 
fleet. 
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Creating a file 

The “New” command in the File Menu provides the option 
to create any of the file types enumerated above. 
Depending on the type of file to be created, the data entry 
user interfaces will be automatically loaded. The “Open” 
command in the File Menu allows the loading of any of the 
file types and users have the option to load the same file as 
many times as desired, depending on the capacity of the 
computing environment. Loaded data can be edited and 
saved. Care should be taken when the same file is loaded 
more than once not to overwrite previously edited records. 

All inputs in the memory of the computer can be uniquely 
catalogued on the designated storage medium by providing 
a filename. In the Microsoft Windows environment, the 
number of characters making up a filename is not limited to 
only eight characters, as in the case of IBM DOS files, and 
users are encouraged to use more descriptive filenames to 
facilitate file identification and retrieval. 

Though the choice of filenames is left to the user, a 
coherent system for naming files is still advisable. One 
practical approach for naming files (now commonly 
applied) is through a fixed format, e.g., GSSSAR. G stands 
for the first letter of the fish's family name, SSS for a 
combination of the first letter of the genus and the first two 
letters of the species name, and AR for the area or stock 
from which this file was taken. Using this format for a data 
set on Upeneus moluccensis, the file may be named 
MUMORG where M refers to the family name (Mullidae), 
U from Upeneus, MO from moluccensis and RG for Ragay 
Gulf, Philippines. 

This format can be further expanded to allow for further 
identification. For example, if the length-frequency data in 
MUMORG, a file originally with 1 cm class intervals, were 
regrouped to 2 cm, it would be appropriate to keep the 
original file and to use another filename for the new file. 
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The new filename may contain the original name 
MUMORG plus a suffix, say 2, referring to the larger class 
interval. Thus, the new filename would be MUMORG2. 


WARNING! 

New data entered or edited are not automatically saved on 
disk. The "Save" and “Save As...” option under the FILE 
menu of FiSAT II must be used to save data that have been 
entered. It is good practice to save a file right after it has 
been created or edited to avoid losses, e.g., due to power 
failures. 


The DATA subdirectory of FiSAT is the default address 
when saving and retrieving files. However, FiSAT II 
provides the option to save data to a different address. 
However, as with filenames, it is advisable to also label 
data folders properly for efficient retrieval of records. 


WARNING! 

Edited data are not automatically saved to disk. 
THEREFORE, IT IS ADVISABLE TO FIRST SAVE THE 
DATA AFTER EDITING BEFORE ANY ANALYSIS IS 
PERFORMED. 


Printing a file 

The "Print..." option of the File Menu allows the user to 
print the loaded file to a printer or to a file. The “Print...” 
command is only applicable to data files and not to results 
of data analysis. A separate “Print” command button is 
available to users if the results are to be printed or graphs 
saved to disks. 

The output device can be reconfigured from the “Print 
Setup” command of the File Menu. The printer re- 
configuration option can also be accessed from the Printer 
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Control Panel of Microsoft Windows or from a command 
button in the “Print” dialog box. 


Sample manipulation 

As defined in earlier sections, file manipulation pertains to 
file operations which alter the original content of the file. 
The changes which result from these manipulations are not 
automatically stored to disk. It is therefore necessary to 
save the altered file if the new data need to be accessed for 
later use. 


WARNING! 

If the manipulation of the data generates values greater than 
the upper limit or smaller than the lower limit set for a 
field, FiSAT II will display a warning message. Also, the 
high value(s) will be replaced by the upper limit (and 
conversely for low values). 


A number of routines are presented below, with 
information arranged as follows: 
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Required file 

Defining the type of file required; 

Input parameters 

Listing the parameters required 
from the user; 

Function(s) 

Defining the mathematical model 
or function used by the routine and 
identifying its variables; 

User interface 

Providing notes on how to use the 
routine. 

Remarks 

Providing some comments and 
remarks to better describe the 
routine; 

Readings 

Providing a few suggested readings 
on the principles behind the 
routine (arranged in alphabetical 
order ) other than the required 
readings, Gayanilo, Sparre and 
Pauly (1996), Gayanilo and Pauly 
(1998) and Sparre and Venema 
(1992). 


Merging two files 

Similar files may be merged, resulting in consolidated files. 
FiSAT II will not prevent users from merging several files 
with different labels. In such cases, however, the user 
should operate with care (as in any other file manipulating 
routines). 

Required file Length-frequency or weight-frequency 

data file with constant class size. 

Input parameter Source files to merge (file already 
loaded into memory as default input) 
and destination file. 
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Function None 

User interface This routine (Fig. 3.16) will require 

only the identification of two length 
frequency files with the same class sizes 
and filename where the merged file will 
be saved. The file loaded is by default, 
one of the two files to merge. However, 
the inputs for the source files can be 
altered. 


*J 


Metge tht He 
vdh fie 


Saw to. 


|C AProgram Fies\FiSAT HNDataVSampleVq 
jCAFiogram F4«\F(SAT H\0a»a\Sample1 Iq 


|C \P 109 am nkn\FiSAT ll\Data\S*npte2 Iq 
proceed | £ancd 


_J 




Fig. 3.16. User interface in merging two length frequency files. 


Remarks The same file may be merged to itself. 


Correcting using probabilities of capture 

Length or weight frequencies biased by gear selection can 
be (partly) corrected using this routine, which divides 
frequencies by their corresponding probabilities of capture. 

FiSAT II provides two options (Fig. 3.17), (1) using a 
probability file or (2) using selection parameters (L 50 , L 75 , 
R5O’ an d ^75)- 

Required file Length- or weight- frequency data file 

with constant class size 
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Input parameters Option 1: Filename of the probability 
of capture file 


Option 2: L 50 , L 75 , R 50 , and R 75 


In both cases, the destination filename will be required. 


j Correct L/F for Probabilities.. 


Coned this We: |c.\Ptcgram FilesVfiSAT ll\Data\Satnpte Iq 

Options 

6 Using a piobaMily Me 
PrcbaStWy file f~ 

Usng selection parameters 

i c ' | n75 | roooo 

-• | bso | wtnn 

Save corrected fiie to: ] 

Proceed | Cancel 

Fig. 3.17. User interface for correcting length frequencies for probability 
of capture. 

Functions The new frequency (N new ) for the class 

group L, is computed from 

Nnew — N 0 |,j / P L 

where P L is the probability of capture. 
If Pl is to be computed (i.e. Option 2), 
this is done using: 

P,=l /(I +exp(S.-S 2 *L))v.. 

( I + exp (S 3 -S 4 *L)) 

where 



S, =((L 50 -ln(3))/(L 75 -L 50 )), 
S 2 : Si/L 50 , 
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5 3 “ (( L 50* ln (3)V( R 75‘ R 50))’ 

5 4 = S3/R-50 

User interface The option to use in correcting the 
selected length frequency file for errors 
associated with probability of capture is 
presented in radio buttons. To select an 
option, click the button to enable the 
other input fields. Note that the default 
values for R50 and R75 are huge values 
indicating there is no deselection. 
However, these values can be altered. 

If Option 1 is used and the class size of 
the file which contains the selection 
data does not coincide with the length- 
or weight-frequency data file loaded 
onto memory, FiSAT II will cancel the 
operation and prompt for another 
filename. 

Remarks Setting the right-hand selection (or 

"deselection") to very high values 
would have the effect of ignoring 
deselection completely for Option 2. 


Adjusting class intervals 

The routine to adjust class sizes can be used to (i) allow 
two files with different class sizes to be merged, (ii) reduce 
the number of classes to fit into the limits set by FiSAT 11 
or (iii) smooth out irregularities, and hence improve the 
interpretability of a data set. In FiSAT II, this routine 
regroups the frequencies by "slicing" the original file. 

Required file Length- or weight-frequency data file 

with constant class size 

Input parameter New class size 
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Functions 


User interface 


The frequencies of the new classes have 
the value: 


N 


new 


= Ni+ N 2 


where 


N, 


-Nj 


f x i4 , 


-Y, 

-Xj 


N 2 -N, tl 


Y — Y 
H+l A i+1 

k ^i+2 “ ^i+1 , 


where 


Nj - frequency of the first class, 

N i+ i - frequency of the following 
class, 

Xj - lower limit of the first class, 
X i+1 - upper limit of the first class and 
lower limit of the following 
class, i+1, 

Xj +2 - upper limit of class i+1, 

Yj - lower limit of the new class, 
and 

Yj +] - upper limit of the new class. 

As a default, the active length 
frequency file will be the default 
source file. A new class size and the 
destination file will be required by the 
routine. 
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xj 


Swrcthle |C \Progr«m Ftos\F£AT IINOataVSampIt.Vq 

Fie P««meta* 

Species name |UpmM rnoluccenw 
Other He ndentfer: |RAGAVGJLF I960 

Class ice 

Old class toe I 1 00 New class tae \~ 


Dest nation tie: 


J 


froceed £kwe 


Fig. 3.18. User interface for adjusting the constant class size of a length 
frequency file in FiSAT II. 

Remarks Adjusting class sizes is another 

approach for smoothing the data files. It 
is advisable to always use the original 
file when adjusting class intervals and 
not to use files previously manipulated. 


Smoothing data files 

This routine is meant to smooth out irregularities (Fig. 
3.19). Smoothing can be done by: (1) a running average 
over 3 classes or (2) a running average over 5 classes. 


Required file Length- or weight-frequency data file 

with constant class size 


Input parameter None 

Functions Option 1 : 

Nj.ncw = (Nj-i + Nj.old + Nj+j )/3 

Option 2: 

N i,new = ( N i-2 + N i-I + N i,old + N i+1 + 

N i+2 )/5 
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■JM W HW W I MWWB 


Source lile: |CAProgram FdesVFiSAT ll\Dala\Sample la 


XJ 

J 


File Pet emeteis 


Species name' |upeneus moluccensis 


Olhei file indenlifier |RAGAV GULF I960 


Smoolhir»g option 

(* Running average over 3 dawet Running average over 5 classes 


Destination file 


J 


Proceed 


Close 


Fig. 3.19. User interlace in smoothing length frequencies using running 
averages. 


User interface 


Remarks 


The default source file is the active 
length frequency file. Simply click the 
desired option, and provide the inputs 
for the desired destination file to 
complete the inputs. 

For both options, pseudo classes with 
frequencies of zero (0) are added below 
and above the extreme sizes, i.e. the 
smallest and largest classes before 
computing the new frequencies. 
FiSAT II will not prevent users from re- 
applying the smoothing function. 
However, this is not recommended, as 
repeated smoothing leads to loss of 
information. 


Readings Laurec and Mesnil (1987) 

Pooling samples 

This routine (Fig. 3.20) allows a user to pool samples from 
several sampling dates. The weighting options (by the 
square root of sample total or % of sample total) that must 
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be considered before actual pooling are presented further 

below. 

Required file Length- or weight-frequency data file 

with constant class size 

Input parameters Identification of samples to be pooled. 

Function N new = N sam p] e j + N sam p| e 2 + ... 

+ ^sample n 

User interface As with previous routines, the default 
source file is the active length 
frequency file. Using the mouse, select 
the samples to pool, provide the new 
sampling date (default is the average of 
the dates of the samples selected). 
Provide the address of the destination 
file and click the command button 
“Pool Selected Samples”. Continue the 
process until all samples to pool has 
been completed. Click the command 
button “Save File” to save the results. 


Source 

Filename 

Specie* 


|C SPro <pan> FtetVFISAT IIVOa»e\Sanv*e lq 


Otoe* •dentate [RAGAY GULF I960 


Dettanabon hie |C kPtoqram FtetVFISAT IIVBataVPooto dS aynple* lq 
Select Samples to Pool 



Enter New Samping Date 

1 22/OS/I 900 


Sample 5 (27/09/1 9801 
I Sample S (20/11/19801 
I S«rpte 7(11/12/1900) 


—3 

Pool Selected Sampler | 


Save Fie and E«r 


Cancel 


j*l 




Fig. 3.20. User interface in pooling length frequency files. 


Remarks It will be recalled that for length or 

weight frequencies stored through 
FiSAT II, a header may be attached to 


Copyrighted material 



43 


each sample (e.g., sex, type of gear, 
depth, etc.). When pooling samples, 
these values (or entries) will be retained 
if and only if they are the same. Entries 
for total weight of the sample, total 
weight of the species and total weight of 
the catch will be added. 


Converting files 

Certain data sets need converting before they can be used. 
FiSAT II (Fig. 3.21) provides two options using cither: 
L ntw = d + e • L old ' (or length-to-length) or L ncw = d • W c 
(or weight-to-length) functions. 



Fig. 3.21. User interface in converting grouped frequencies in FiSAT II. 

Required file Length-frequency data file with 

constant class size 
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Input parameters Depending on the option selected the 
parameters for d, e and f (for 
Option 1 only). 

Function Option 1 : L nC w = d + e • L 0 i/ 

Option 2: L new = d • W c 

User interface Enter the address of the destination 

file, select the conversion function to 
use using the mouse and enter the 
parameters before proceeding. 
Option 1 is used for converting one 
length measurement into another, 
e.g., total length into fork length or 
carapace length to total length when 
the variable “f” is set to 1 . 

Remarks Option 1 can also be used for 

conversions based on a 
proportionality (e.g., inches to cm), 
by setting d = 0, and e = conversion 
factor. Also, a constant can be added, 
by setting e = 1 and d = constant. 

Option 2 can be used to convert 
weight-frequency data to length- 
frequency data given the parameters 
(a,b) of a length-weight relationship 
by setting d=(l/a) (|/b) , and e=l/b. 

Raising samples 

For some of the fish stock assessment methodologies 
included in FiSAT II, length or weight frequencies may be 
required to represent certain quantities, for example, total 
catch. The following (Fig. 3.22) are options for raising 
frequencies: 
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Fig. 3.22. User interface in raising length frequency files. 


Using catch per unit of effort (C/f) 

When the total catch and a measure of effort is available, 
the frequencies can be raised to represent the catch per unit 
of effort (C/f). 

Required file Length- or weight-frequency data file 

with constant class size 

Input parameters Catch per unit of effort (C/f) per sample 

Function N new = N 0 i(j • (C/f) 

User interface Select the samples to raise and click on 
this option to select. 

FiSAT II will use the C/f data provided 
on the sample sub-header as default 


Copyrighted material 


46 


values; in the absence thereof, FiSAT II 
will cease processing. 

Remarks In a multifleet fishery, data may need to 

be raised using this option to obtain a 
better representation of total catches. 


Using catch/sample weight 


When in addition to the total sample weight (WS), the 
weight of the total catch (WC) is available, this option will 
allow frequencies to be raised to better represent the total 
catch. 


Required file Length- or weight- frequency data file 

with constant class size. 


Input parameters Total weight of the catch (WC) and 
total weight of the sample (WS), or the 
'a' and 'b' coeficient of the length-weight 
relationship to estimate total sample 
weight. 

Function N new = N 0 | ( j*(WC/WS). 

User interface Select the samples to raise and click on 
this option to select. 

FiSAT II will use the WC and WS data 
provided on the sample sub-header as 
default values. In the absence of 
required parameters, the routine will 
prompt an error and will discontinue 
processing. 

Remarks This option may be necessary to raise 

data to represent the total catch required 
for the length-structured VPA and 
related routines. 
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Using a constant 

This option may be used when independent estimates of an 
appropriate raising factor (RF) are available, or to change 
units of measurements, e.g., inches to cm. 


Required file Length- or weight-frequency data file 

with constant class size 

Input parameters The factor (0 to 999,999) to raise the 
frequencies. 


Function N new = N 0 | d • RF 

where RF is a constant 


User interface Select the samples to raise and click on 
this option to select. 

When this option is selected, FiSAT II 
will prompt for input of the constant. A 
default value of 1 (one) is provided. 

Remarks None 


Percent of sample total 

This option allows the frequencies of samples to be re- 
expressed as % frequencies; this option can be used to give 
equal weight to all samples of a file; preferably it should be 
applied to samples having similar sizes to start with. 

Required file Length- or weight-frequency data file 

with constant class size. 

Input parameters None. 

Function The new frequencies are defined by 

^new = [N 0 i d /Z(N old )] • 100 
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where: 

N new = new frequency, 

N 0 ld = original frequency, 

2(N 0 id) = tota l frequency in the sample. 

User interface Select the samples to raise and click on 
this option to select. 

Remarks This type of data transformation is also 

available in some of the routines that 
require data to be pooled (e.g., catch 
curve analysis); however, these routines 
do not allow saving of the data thus 
transformed. 


Square root of sample total 


This option allows the frequencies of a sample to be re- 
expressed such that their sum becomes proportional to the 
square root of the total sample size; it should be considered 
for samples with initially very dissimilar sizes. 


Required file Length- or weight-frequency data file 

with constant class size. 


Input parameters None 

Function The new frequencies are defined by 

N„cw = [N'old /2(N 0 ,d) ,/2 ] • 100 

where: 


N ncw = new frequency count, 

N 0 )d = original frequency, 

ZN 0 | d = total frequency of the 'old' 
(original) sample. 
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User interface Select the samples to raise and click on 
this option to select. 

Remarks As with the previous option, this type of 

data transformation is also available 
(but without saving) on some of the 
routines that require pooling of samples. 

Temporal weighting 

Temporal weighting will raise each sample by a factor 
proportional to the time separating it from the nearest other 
sample, which gives more emphasis to samples 
representing longer periods and less to samples collected at 
close intervals. 

Required file Length- or weight-frequency data file 

with constant class size 

Input parameters None 

Function N new = [N old • D'fi 

where: 

D is the number of days separating the 
sample from its nearest neighbour. 

User interface (see routines with this prompt: 

e.g. length-converted catch curve). 

Remarks FiSAT II will not weigh samples by 

time if D exceeds 3 months. 
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Importing old FiSAT data 


This routine will allow users of the old version of FiSAT to 
readily migrate to the updated FiSAT II version. 


Required file Length -frequency data file saved in the 

old version of FiSAT 


Input parameters Source and destination files 
Function None 


User interface Enter the source file (i.e. old FiSAT 
data) and the destination file (Fig. 3.23). 

Remarks The contents of the sub-header will not 

be imported other than the data of 
sampling. 



Fig. 3.23. User interface when importing old FiSAT data files to the new format. 

Exporting FiSAT II length frequencies 

In many cases, the FiSAT II file may have to exported to 
other forms for report purposes or to allow data importation 
to other software packages (e.g. Microsoft Excel). 

Required file Length- frequency data file with 

constant class size. 
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Input parameters None 
Function None 

User interface The routine will access the standard 

dialog box (Fig. 3.24) for exporting a 
file. A tab-delimited file will be 
generated that is standard to most 
commercially available electronic 
spreadsheets and word processors. 


rxport... 


Save irt | O FiSAT II T] f B 0’ 



Fig. 3.24. Dialog box used when identifying the destination file of the data to be exported. 


Remarks The sub-headers will not be exported 

other than the date of sampling. 
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Chapter 4. ASSESS menu 


What you will learn from this chapter 

This chapter presents routines for analysis of the various 
types of data presented in the previous chapter. However, 
these are presented in a summarized form, and we remind 
FiSAT users to read the FiSAT reference manual (Gayanilo 
and Pauly, 1998), and other documents explaining the 
background of these routines. 

The format in module description 

The Assess menu is divided into two parts, (i) the routines 
to estimate parameters (e.g., L x , K, C, M, Z, etc.) and (ii) 
the routines to predict yield or stock-related attributes given 
certain fishing scenarios. In the following, the modules are 
presented as given in Chapter 3, with one description 
added: 

Output(s) This section defines the output of the mode I 
or methodology being described. Examples 
will be provided when necessary. 

Parameter estimation 

The step that should immediately follow data entry and 
manipulation is the estimation of population parameters. 
These are required inputs to the second set of analytical 
routines that deals with predictions. 

Direct fit of length-frequency data 

The set of routines classified under this heading are those 
which estimate growth parameters directly from the length 
composition of the stock, without previously translating the 
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length scale into an age scale, as an alternative to the MPA 
approach described below. 

ELEFAN I 


ELEFAN I is a routine that can be used to identify the 
(seasonally oscillating) growth curve that "best" fits a set of 
length-frequency data, using the value of R n as a criterion. 
FiSAT II provides three options to the user to identify that 
"best" growth curve: (1) curve fitting by eye (plotting of 
the histogram or restructured data may also be accessed 
from the Support Menu): (2) scan of K-values (Fig. 4.1), 
and (3) response surface analysis (Fig. 4.2). The plotting 
and curve fitting by eye is described in detail in the Support 
Menu below. 

Note that in ELEFAN I, the parameter t 0 is replaced by the 
coordinates of a point (any point actually) through which 
the curve must pass, and whose coordinates consist of SS (a 
starting sample) and of SL (a starting length). 

Required file Time scries of length-frequency data 

with constant class size. 

Input parameters Option I : Scanning of K-values 

Fixed value of L„ C and WP and, as an 
option, a fixed starting point. 

Option 2: Response surface analysis 

Range of values for two parameters 
(e.g., L* and K), fixed value for the two 
others (e.g., C and WP), and, as an 
option, a fixed starting point. 

Functions In ELEFAN I, data are reconstructed to 

generate "peaks" and "troughs", and the 
goodness of fit index (R n ) is defined by 
R n = IQESP/ASP/io 
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Outputs 


User interface 


where the ASP ("Available Sum of 
Peaks") is computed by adding the 'best' 
values of the available 'peaks' and the 
ESP ("Explained Sum of Peaks") is 
computed by summing all the peaks and 
troughs "hit" by a growth curve of the 
form, 

L t = L„ ( 1 - EXP(-K(t - 1 0 ) + S K + S t0 ) 


where 

S ts = (CK/2jt) • sin (2n:(t-ts)), 

S t0 = (CK/2:t) • sin (2ji(t 0 -ts)), and 
L t is the length at time t. 

Option 1: Scanning of K- values 

This plots R n values for a range of K. 
values (0. 1 0 to 1 0) on a log-scale. 

We recommend use of this plot for aU 
growth analyses, if only to assess how 
reliable an estimate of K is. 

Option 2: Response surface analysis 

This outputs an 11 by 11 matrix 
showing R n values and in which the 10 
best values are highlighted, thus 
enabling selection of the "best" 
combination of growth parameters. 

ELEFAN I contains three tabs. The first 
tab is only to identify the length 
frequency file to use. The tab, “K Scan” 
(Fig. 4.1), would allow the user to view 
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the behaviour of the scores given 
estimates of L*. 

The starting point can be set to variable 
or fixed starting point (see red arrow in 
Fig. 4.1). If starting point is set to a 
fixed value, the starting sample and 
starting length can be selected from a 
drop-down list. 


JjELFFAN I [C:\Pro<jram Ne*\RSAT UyData Sample IfqJ... 


NonPaiametnc Scoring of VBGF Fit Using ELEFAN I 

general & Scan j Response Surface] Automatic Search] 


*1 


1 60 



2 30 2.60 

Growth performance Index (•’) 


Enter Parameter* foe K-scan 
Loo: 21 

C: 000 

WP 000 

W Fix the rtartmg port 

Starting sample fl 
Starting length. 1 9.00 H 


3 30 3 60 


Parameter* at Maxrmu 

m 

Loo 

21 00 

K: 

0 760 

C: 

0 000 

WP 

0000 

Starting sample 

1 

Starting length 

900 

Score 

0415 



/■ 

V. 

P Ahrend « 


viLl 

xl 


Plot VBGF 


Curve | Compute | grint | £to*e 


Fig. 4.1. User interface of ELEFAN I. K-scan routine. The red line in the plot is the 
actual score and the thick blue line is the overlaid trend line. 


The behaviour of the scores may be so 
erratic that assessment of the possible 
solution is difficult. A trend line (see 
red circle in Fig. 4.1) may be applied. 
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The scroll bar is used to change the 
resolution of the trend line.The third tab 
is to execute a response surface analysis 
(Fig. 4.2). The starting point is a 
requirement of this subroutine. Note 
that the value of the L* cannot be less 
than the starting length. 


•■FI FFAN I [C:\Program File*\RSAT II\Data\ 




N on-Parametric Scoring of VBGF Fit Using ELEFAN I 


*1 


general) £ Scan ReJPonie Surface j Automate S patch] 
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Parameters 

From 

To 

parameters T o make a parameter constant. 




enter the same value for lower and upper 




kmit 
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1 
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000 

- 0.00 
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WP 

000 

0.00 


Scores ELEFAN I Method 
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20.20 
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0.185 
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0.74 

0.213 

0.227 

0.211 

0.211 
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0.227 
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0 255 
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0.82 

0.261 
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0.310 
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0.310 
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0227 

0 86 

0.272 
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Plot VBGF Curve | Compute | Erint | Close 


Fig. 4.2. User interface of ELEFAN I, Response surface analysis. Note the change of 
colour indicating improvement on the score. 


Remarks It is imperative that users of ELEFAN I 

read detailed accounts of this method. 

Readings Pauly ( 1 982) 

Pauly and David (1981) 

Pauly and Morgan (1987) 


Copyrighted material 


57 



Conceptually, this approach is similar to ELEFAN 1 in that 
it is designed to maximize a non-parametric scoring 
function. Two options for identifying optional values of L, 
and K are available: (1) response surface analysis and (2) 
scan of K-values, both of which are very similar, in display 
and operations, to those of the ELEFAN I routine (see 
above). 


Required file Time series of length-frequency data 

with constant class size. 


Input parameters Option 1 : Response surface analysis 
Range of values for L, and K. 
Option 2: K-scanning 
Fixed value of L*. 


Functions The score (S) for Shepherd’s method is 

defined by: 

S = (s A 2 + sb*)* 

where s A and s B are the goodness-of-fit 
scores (s tz ) obtained with the origin of 
the VBGF in calendar time (t^ set to 0 
and 0.25, respectively. s tz is defined by: 

su-^t.-vn; 

t 
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Outputs 


User interface 


where 

Nj = frequency for length group i, 
Tj = D • cos 2 jt (t-tj), 

D = (sirnr (At)/jr(At)), 
t = At/2, 

^ — Viax " tmin> 

tj — t z - (l/K)-ln(l-(Lj/L*)), and 

t z = (l/2jt) ■ tan'*(s B /s A ). 


Option 1: This outputs an 11 by 11 
matrix of S values (with S max 
standardized to 1) and in which the 10 
best scores are highlighted, thus 
enabling selection of the "best" 
combination of L, and K; 

Option 2: This generates a plot of S 
values (with S max standardized to 1 ) for 
a range of K. values (0. 1 to 1 0 year 1 ) on 
a log scale, thus enabling the 
identification of the best value of K for 
a given value of L». 

The user interface of this routine is very 
similar to that in ELEFAN I. It also 
contains three tabs where the first tab is 
used only to identify the file to be 
analysed. The second tab (Fig. 4.3) is to 
examine the S values given a fixed L„ 
and a fixed range of K values. 
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Fig. 4.3. User interface of Shepherd's method. The trend line indicates a K value 
in the same range as that generated by ELEFAN I in Fig. 4.1. 


The third tab (Fig. 4.4) is the response 
surface analysis given a range of L* and 
K values. The main difference with 
ELEFAN I is that, in this approach, the 
starting point and seasonality 
parameters (C and WP) are not 
identified. 
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Fig. 4.4. Results of the response surface analysis where red (i.e. darker coloured 
cells) indicates improvement on the score. 


Remarks As for the ELEFAN I routine, this 

method is best applied when L* has 
been estimated using another method 
(e.g., the Powell-Wetherall plot, see 
below). Note however, that this method 
differs from ELEFAN I in that it cannot 
deal with seasonal growth oscillations. 
Also, note that in this implementation, 
the score function is standardized to 1, 
i.e. S max is made equal to unity. 
Further, the parameter "t z " (which is not 
similar to t 0 ), is here replaced by a 
"starting point" (as for ELEFAN I, see 
above), and hence the ELEFAN I output 
routine can be used to display growth 
curves whose parameters have been 
estimated using Shepherd's method. 

Readings Shepherd ( 1 987) and Isaac ( 1 990) 


Copyrighted material 



61 


Powell-Wetherall plot 

This method allows estimation of L* and Z/K from a 
sample representing a steady-state population, as can be 
approximated by pooling a time series of length-frequency 
data. 

Required file Length-frequency data file with 

constant class size 


Input parameter Graphical identification of smallest 
length fully recruited by the gear (L\ or 
cut-off length): 


Functions 


(L-L') = a + b • L’ 


where 


l f L. +L' 

\ 1 + (Z/K.) 

from which 

L» = -alb, 

Z/K = - (l+b)/b. 

Outputs L, and Z/K estimates and regression 

parameters. 

User interface In this routine, as with other routines in 
FiSAT II, options are provided to select 
and deselect samples to be included in 
the analysis. To deselect a sample, 
uncheck the sample represented by the 
sampling date. Three options are 
provided to do a temporary data 
transformation (see Chapter 3 for 
details). Use the radio buttons (Fig. 4.5) 
to select the option. In addition to the 
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available temporary data transforma- 
tion, temporal weighting may also be 
applied to the selected samples by 
checking the appropriate box (see 
Chapter 3 for details of this option). 



Fig. 4.5. File identification tab and options to temporarily transform data before 
data analysis. These options are available in other routines of FiSAT II that sums 
a time series of length frequencies to approximate a steady-state condition. 

Once the initial inputs are satisfied, the 
data points may be plotted (Fig. 4.6). In 
FiSAT II, the identification of the cut- 
off length (L') is facilitated by an 
adjacent graph (pseudo-catch curve; 
Fig. 4.7), showing one's L/F data, and 
from which L' is often identifiable. The 
cut-off length may be identified using a 
mouse pointing device directly from the 
Powell-Wetherall plot or from the 
pseudo-catch curve to estimate L x and 
Z/K or from the pseudo-catch curve. 
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A command button, “Open pseudo- 
catch curve” is available in Fig. 4.6, if 
the pseudo-catch curve (to help identify 
the cut-off length) was previously 
closed. 



Fig. 4.6. The Powell-Wetherall's plot showing the regression line with the Y- 
intercept as the estimate of L». 
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B *J 

NOTE T hi* pto* n a pseudocatcfi curve to hdp idenffy the age 
group fidy recruted lo the fishery and we (age) not affected by gear 
setecbon Identfy the group v*th a mouse dick 


PSEUDO-CATCH CURVE 



Close fat Wndow 


Fig. 4.7. Pseudo-catch curve to help identity cut- 
off length in Powell-Wetherall's plots. Selecting the 
cut-off length can also be done through this plot. 


Remarks This approach should be used only 

when sufficient samples are available. 

Readings Pauly ( 1 986) 

Powell (1979) 

Wetherall (1986) 

Analysis of length-at-age data 

This routine allows non-linear estimation of growth 

parameters from length-at-age data. 

Required file Length-at-age data, consisting of at 

least four age-length data pairs. 
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Input parameters Option 1: Non-seasonal growth 

The minimum and maximum 

constraints and starting values for the 
growth parameters L x , K, and t 0 . 

Option 2: Seasonal growth 

The minimum and maximum 

constraints, and starting values for the 
growth parameters L x , K, t 0 , C and WP 
(where WP = tg + 0.5). 

Functions Find the global minimum of the 

following functions using Marquardt's 
algorithm: 

Option 1 

SSE - ^ i. " L *[' - EXP(-K(t, - 1 0 ))} 2 


Option 2 

SSE - J i, - L* D - EXP(-K(tj - 1 0 ) + S B + S l0 >}’ 

I 

where 

S ts = (C‘K/2 ji) • sin (2ji(t-t s )), 

S t0 = (C-K/2 jt) • sin (2;r(t 0 -t s )). 

Outputs L x , K, t 0 , and their standard errors for 

Option 1 and in addition, C and WP for 
Option 2. For either option, the 
coefficient of determination (r 2 and 
adjusted r 2 ), phi prime (o'), and its 
standard error. 
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User interface The user interface of this routine 
contains only two tabs (Fig. 4.8 and 
Fig. 4.9). The first tab (Fig. 4.8) 
provides the option to enter the range of 
values identified by the minimum and 
maximum constraints, and the starting 
estimates for each parameter. After the 
iteration process, the results are also 
displayed in this tab. 



Fig. 4.8. The "Parameters" tab of the user interface when analysing length-at-age 
data. This tab identifies the file and range of values to search for the growth 
parameters. The results, including the statistics to measure the fit, are also 
displayed after iteration. 
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The second tab (Fig 4.9) is a plot of the results, i.e. the data 
points overlaid by the computed growth curve. 



Fig. 4.9. Plot of the growth curve based on the computed growth parameters. 


Remarks For rapid convergence, the routine 

should be initiated with a wide range of 
constraints. To fix inputs at certain set 
values, enter minimum constraint = 
maximum constraint. 

Note that up to 200 iterations may be 
performed until the routine converges or 
stops. 

The method for estimating the standard 
error of o' stems from J. Hoenig (pers. 
comm.). 
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Readings Kvalseth ( 1 985) 

Pauly et al. (1992) 

Saila, Recksiek and Prager (1988) 
Soriano and Jarre ( 1 988) 

Modal progression analysis 

Modal progression analysis (MPA) refers to a methodology 
that infers growth from the apparent shift of the modes or 
means in a time series of length-frequency samples. MPA 
involves three stages: (i) decomposition of composite 
distributions into their components to identify means, (ii) 
subjective identification and "linking" of the means 
perceived to belong to the same cohorts and (iii) using the 
growth increments and size-at-age (relative) data resulting 
from the linking to estimate growth parameters. 

In FiSAT II, two methods are provided to decompose 
composite length-frequency distributions^) Bhattacharya's 
method and (ii) NORMSEP. 

Bhattacharya's method 

Required file Length or weight frequencies with 

constant class size. 

Input parameters Visual identification of frequencies 
perceived to belong to one age group, 
using the graphs provided for this. 


Functions In (Nj + |) — In (Nj) = aj + bj-Lj 


where 

Nj and N;+j are successive frequencies 
of the same component of a group of 
fish in a sample (i.e. representing age 
group j) and where Lj is the upper class 
limit of Nj. 
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From this, the mean of the normal 
distribution is 

Lj = -aj/bj 

while its standard deviation (a) is 
oj = (-AL/bj) 1/2 

where AL is the constant class size. 

Also, a separation index (SI) is 
computed. 

SI = ALj/AOj 

where ALj is the difference between two 
successive means, and Aoj is the 
difference between their estimated 
standard deviations. 

The separation of length-frequency 
samples into their components is an 
iterative process in that every identified 
component is subtracted from the 
remainder of the sample using the 
Gaussian function. 



where Nj + refers to the previous set of 
frequencies and N 2 j+ is the new set of 
frequencies, less the component 
identified so far. Also, the user can 
backtrack and redo any step of the 
analysis without losing intermediate 
results. 


Copyrighted material 



70 


Outputs 


User interface 


Mean lengths, population sizes (in 
numbers), standard deviations and 
separation index (SI) for the age groups 
identified. 

These results are automatically saved to 
disk as a "mean and standard deviation" 
file, using the same filename as the 
input. 

The module contains three tabs 
(Fig. 4.10). As with previous routines, 
the first tab is used only to identify the 
file to analyse. The “Plots” tab shows 
the log-plots of the slopes between 
successive frequencies and the resulting 
distributions. 

The process starts with the selection of 
the sample to analyse (sec red circle in 
Fig. 4.10) and identification of the first 
and last data points that identify a group 
using the mouse. 
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Fig. 4.10. User interface for Bhattacharya's method in FiSAT II. 


Clicking the “Redo” command button 
will re-initialize the display and the user 
will be provided with the option to 
reselect data points. Clicking the 
“Continue” command button indicates 
acceptance of the results and the group 
will be eliminated from the composite 
distribution to emphasize the next 
possible group (if any). This process is 
continued until all groups are identified. 
The result can be saved (Fig. 4. 1 1 ) by 
clicking the “Save MSD” command 
button that replaces the “Continue” 
command button after group 
identification. The results will be saved 
as a "mean and standard deviation file" 
that can later be linked to form growth 
increments or saved as length-at-age 
(relative) data. 
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Fig. 4.11. "Save MSD" command button replaces the "Continue" command 
button in Fig. 4.10 every after successful identification of an age group. 


A summary of the results, including the 
SI values, is presented in the third tab 
(Fig. 4.10). 

Remarks Up to 10 age groups can be identified 

per sampling period, but separation is 
generally unreliable when the SI value 
is below 2. Also, note that the standard 
deviations usually increase with lengths. 

Reading Bhattacharya (1967) 

NQRMSEP 


This method applies the maximum likelihood concept to 
SEParation of the NORMally distributed components of 
size-frequency samples. 

Required file Length or weight frequencies with 

constant class size. 

Input parameters The expected number of age groups 
(maximum of ten) and their expected 
mean lengths. [The results from 
Bhattacharya's routine (see above) are 
used as initial inputs when available.] 
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Functions 


Outputs 


User interface 


(Refer to suggested readings for the 
maximum likelihood method) 

Mean lengths, population sizes (in 
numbers), standard deviations and 
separation index (SI) for the groups 
identified. 

These results may also be saved to disk 
as a "mean and standard deviation file" 
using the same filename, and 
overwriting results from previous 
analyses. 

The user interface of this routine is very 
similar to Bhattacharya’s routine. 
However, it contains only two tabs (see 
Fig. 4.12). The first tab identifies the 
file to analyse and the second tab is for 
group identification and a summary of 
the results. 

Users first select the sample to analyse 
then approximate the mean length of 
each group (normally, the results of 
Bhattacharya’s routines are used as 
defaults). Clicking the “Compute” 
command button iterates for the best fit 
using the approximate means as a 
starting point. The computed mean 
lengths can be saved by clicking the 
“Save MSD” command button as in 
Bhattacharya’s method. 
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Fig. 4.12. NORMSEP user interface in decomposing composite distributions. 

Remarks It is advisable to use Bhattacharya's 

method to obtain initial guesses and to 
refine the results using NORMSEP. 

NORMSEP, once initiated, will run 
until it reaches an optimum, or 200 
iterations. 

Readings Abrahamson (1971) 

Hasselblad ( 1 966) 

Pauly and Caddy (1985) 


Linking of means 

Required file Mean lengths and standard deviations 

file. 

Input parameters Links between means thought to 
represent the same cohort. 
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Function None: linking of means results from 

user's subjective interpretation of a 
graph. 


Outputs Either (1) a number of links between 

means (i.e. growth increment data), for 
analysis with a Gulland and Holt plot, 
then refined using another method, or 
(2) series of linked means, stored as a 
length-at-age (relative), for analysis 
with the appropriate routine. 

User interface The module requires a mean and 
standard deviation file (*.MSD) that 
can be created from the results of 
Bhattacharya’s analyses or from 
NORMSEP. The file can also be 
created manually (see Chapter 2). The 
points can be linked using the mouse 
(Fig. 4.13). 



Fig. 4.13. Linking of mean lengths in FiSAT II. 
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Clicking the subsequent data points 
automatically links the points. To break 
the link, i.e. creating a new cohort, 
double-click in any part of the plot. 
Clicking on the “Save Links” command 
button will save the links as growth 
increments or as length-at-age (relative) 
data. To redo the links, click the 
“Redo” command button. 

Remarks Linking of the mean lengths without 

any prior knowledge of the growth 
patterns of a given species is usually 
extremely difficult, and should not be 
attempted if the mean lengths do not 
display a clear progression through 
time. 

The most likely of two alternative 
hypotheses concerning the distribution 
of lengths about age are: 

(i) constant CV, or 

(ii) CV changing with AL/At and hence 
with length. 

However, for the second hypothesis to 
be likely the analysis should indicate a 
significant correlation; also, the CV of 
L* must be positive. 

Reading Sainsbury ( 1 980) 

Analysis of growth increment data 

FiSAT II provides four options to estimate the growth 
parameters of the von Bertalanffy Growth function (VBGF) 
from growth increment data, which may originate from 
linking of mean lengths (see above), or from 
tagging/recapture experiments. 
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Gulland and Holt plot 

This method allows preliminary estimation of growth 
parameters from growth increments, based on the fact that 
under the VBGF, growth rate declines linearly with length, 
reaching zero at L,. Also, the residuals of the plot are used 
for inferences on seasonality of growth. 

Required file Growth increment data file with at least 

3 pairs of points. 

Input parameters Data points to be included in the 
regression analysis; also independent 
estimate of L„ if such value is to be 
used (in a "forced" Gulland and Holt 
plot). 

Functions AL/At = a + b • L, 

where 

AL - L r - L rn , 

At = t,. - t m , and 
L = (L r +L m )/2. 

where L m is the length at marking 
(initial reading), L r is the length at 
recapture, and t m and t r the 
corresponding dates. 

From these, growth parameters can be 
estimated using 

L x = -a / b and 
K = -b 

Also, the plot of residuals, re-expressed 
in %, and plotted against the midrange 
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Outputs 


User interface 


of the time at large are used for 
preliminary identification of growth 
oscillations, through an iterative routine 
which identifies the time of the year that 
maximizes the difference between the 
means of two sets of residuals, covering 
6 months each. A t-test is used to 
identify significant differences, and a 
preliminary estimation of C is obtained 
from: 

C ^ Ad (n/2) 
where 

Ad relates the difference of the means of 
the two sets of residuals to the 
amplitude of the sinusoid growth 
oscillations. 

Preliminary estimates of K and L x , or of 
K alone if a "forcing" value of L x is 
used; a t-test is used to indicate whether 
the residuals of the plot can be 
separated into two significantly 
different groups pertaining to two 
periods of sixth months, and hence 
suggest seasonal oscillations. If the test 
is positive, preliminary estimates of C 
and WP are outputted. 

The routine contains three tabs (Fig. 
4.14). The parameters (L» and K) are 
automatically computed once a file has 
been identified. 
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Fig. 4.14. User interface for the Gulland and Holt plot. A fixed value of L. can be 
entered (red circle). 


The user also has the option to encode a 
fixed value of L„. In which case, only 
the VBGF growth constant, K, will be 
computed. Growth seasonality will be 
inferred from the residual plot (Fig. 
4.15) presented in the second tab based 
on the computed growth parameters. 
The third tab of the module summarizes 
the statistical results of the regression 
analysis to estimate the growth 
parameters from the Gulland and Holt 
plot. 


Copyrighted material 


I 


80 




* J 


Gulland and Holt Plot 


Header 

Fiename 


|c AProgwn FiesSFiSAT ll\Data\Sample1 gr 


Species name: |La>geRsh 




Other identtfiers |New mafoc exercise (trawler) 
fi&HPIol Residual P w]r egression Analysts | 



NOTE The analysis of the 
residuals ndcates a strong 
osciation o t data points 
rrdcating seasonably on the 
growth of the specimens 


Seasonal growth 
oiciations are stgmhcant 
but small (C<0 101 


Compute with Fixed loo 

Print 

Qose 




Fig. 4.15. Residual plot of the Gulland and Holt plot indicating seasonality 
of the growth. 


Remarks 


Readings 


An option is provided to omit selected 
data points from the analysis, and/or to 
set the value of L„ leading to a "forced 
Gulland and Holt plot". 

Gulland and Holt (1959) 

Pauly (1984a) 


Munro's method 


This method, based on Munro (1982), uses growth 
increment data to estimate L, and K (Option 1), or K 
alone, given L* (Option 2). 
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Required file 
Input parameter 

Function(s) 


Outputs 
User interface 


Growth increment data file 

Fixed value of L* if fixed value option 
is selected 

Minimize the coefficient of variation of 

ratio = [ln(L.-L m )-ln(L x -L r )]/(t r -t m ) 

where: L m is the length at marking 
(initial reading), L r is the length at 
recapture, and t m and t r the 

corresponding dates. Each growth 
increment leads to an estimate of K, 
given L„, and the variance of this 
estimate is computed; in Option 1 the 
value of L x is selected which minimizes 
that variance; under Option 2 only the 
mean of the K estimate and its standard 
error are output, following input of L x . 

Option 1 : Estimate of K and L, 

Option 2: Estimate of K, given L. 

The user interface of this routine 
contains three tabs (Fig. 4.16). The first 
tab presents the results of the analysis 
that follows immediately after a file has 
been identified for analysis. The user 
has the option to select or deselect data 
points for the analysis through the 
equivalent Gulland and Holt plot using 
the mouse. As points are selected or 
deselected, the results of Munro’s plot 
are also re-computed and the graph re- 
plotted. The regression analysis of the 
Gulland and Holt plot is given on the 
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third tab, “Regression Results from 
Gulland & Holt Plot”. 


2SJ 


Header 

Rename 


Munro's Plot 

|C AProgram ResVFiSAT ll\De»a\Sample1 grt 


Species name | Large fish 


Olher identifiers |New mafor exercise (travel) 
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Loo: 

it r 


22 06 
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Fixed Loo 


I - Use fixed Loo 
Enlei Loo I 


Preil 


£lose 


Fig. 4.16. User interface for the routine using Munro's method. This plot is 
recomputed automatically as data points are selected or deselected from the 
equivalent Gulland and Holt plot. 


Remarks The data pairs with lengths greater than 

L, will be deselected automatically as 
these introduce negative values into the 
logarithmic expression of the function 
Note that this method does not account 
for seasonal growth oscillation. 

Readings Munro ( 1 982 ) 

Pauly ( 1 984a) 
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Fabens' method 

Fabens (1965) suggested a method for estimating L, and K 
by predicting length at recapture (L r ‘) based on the current 
parameter selection and the length at marking (L m ). The 
growth parameters are estimated by minimizing the sum of 
squares of errors (SSE), i.e. the squared differences 
between the observed lengths at second reading (L r ) and 
the predicted lengths (L r '). 

Required file Growth increment data file 

Input parameter Fixed value L Jf fixed value option is 
selected 

Functions Minimize the following functions, using 

the Newton-Raphson iteration method: 

SSE-J(L ri -L „■)? 

where the predicted length at recapture 
(L r ') is given by 

L ri' = L mi + (U-L mi )(l-exp(-K.Atj), 
and 

— V " *m 

Outputs Estimates of the L* and K. parameters of 

the VBGF, and their standard errors. 

User interface This module utilizes the same user 
interface as with the Munro’s method. 
The plot (Fig. 4.17) is automatically 
computed after a file has been 
identified. Data points may be selected 
or deselected, using the mouse, from 
the equivalent plot using the Gulland 
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and Holt concept. As points arc selected 
or deselected, the parameters are 
automatically re-computed and results 
re-plotted. The results of the regression 
analysis are given on the third tab. 

Also, users have the option to enter a 
fixed value of L x . 


Fabens' Method— 


Fabens' Method 


2£J 


Header 
Filename 
Species name: 
Other identifiers 


fCAProgram FilesVFrSAT ll\Data\Samplel gin 


| Large fish 
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Compute with Frxed Loo 


Rnnt 


Close 


Fig. 4.17. User interface of Fabens' method in FISAT II. Data points may be 
selected or deselected from the equivalent representation in the Gulland and Holt 
plot. 


Remarks As with Munro’s method (above), 

Fabens ’method does not account for 
seasonal oscillation. 
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Reading Fabens(1965) 


Appeldoom's method 

This method, based on Appeldoom ( 1 987) and Soriano and 

Pauly (1989), allows the use of growth increment data to 

estimate the parameters of a seasonally oscillating version 

of the VBGF. 

Required file Growth increment data file. 

Input parameters The minimum and maximum 
constraints, and starting values for the 
growth parameters L,, K, C and WP. 

Functions Minimize the following function using 

Marquardt's algorithm for a non-linear 
fit: 

SSE - J - (L« - (L« - L, ^XP(- (KAt - S, + S t+At ))f 
where 

S t = (CK/2 ji) • sin (2k (t-ts)), 

St-^ ~ (CK/2 jt) • sin (2k ((t+At)-ts)), 
and, 

ts = 0.5 + WP. 

Outputs Estimates of L„ K, C, and WP and their 

standard errors. 

User interface The user interface of this routine is 
similar to the user interface when 
analysing length-at-age data (Fig. 4.18). 
Users may enter a range of possible 
values and a starting estimate of the 
growth parameters. Default values are 
given based on the data. 
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The resulting growth curve of the 
analysis is plotted in the second tab 
(Fig. 4.19), “Growth Curve”. A 
summary of the statistical parameters 
(partially presented in the second tab) is 
given in the third tab. 


JLAppeldoorn’s Method- 


Header 
Filename: 
Species name 
Other identifiers 


*1 


Appeldoorn's Method 
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Compute 1 1 frint.. | Close 


Fig. 4.18. User interface of Appeldorn's method in estimating growth parameters from 
growth increment data using a non-linear fit approach. 

Remarks For rapid convergence, the routine 

should be initiated with a wide range of 
constraints. To fix inputs at certain set 
values, enter minimum constraint = 
maximum constraint. Note that unlike 
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in the two previous routines, data points 
cannot be selected or deselected from 
this routine. Selection may be done by 
manually editing the source file. 

Readings Appeldoom ( 1 987) 

Soriano and Pauly (1989) 



Fig. 4.19. The resulting growth curve computed using Appeldorn’s model from growth 
increments. 
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Mortality estimation 

Most of the routines to estimate mortality and related 
parameters incorporated in FiSAT II require estimates of 
growth parameters as obtained from the previous set of 
routines. They are grouped in FiSAT II as either (1) 
routines to estimate total mortality (Z) from steady-state 
sample or (2) routines to estimate mortality from natural 
causes (M), here assumed to be constant for all sizes. 

Z from steady-state sample 

The set of routines included here assumes that length- 
frequency samples are available which represent a steady- 
state population (as can be approximated by pooling a long 
series of samples). 

Length-converted catch curves 

FiSAT II provides two options for length-converted catch 
curve analysis: (1) length-converted catch curve for cases 
where growth does not exhibit annual oscillations, i.e. 
C = 0 and (2) length-converted catch curve for cases where 
growlh exhibits annual oscillations, i.e. when C > 0. 

Required file Length-frequency data with constant 

class size 

Input parameters L„ K, and t 0 (t 0 is a facultative entry, 
not required for estimation of Z) (and, 
for Option 2, C and WP) and 
identification of the smallest length 
class that is fully selected and recruited 
(L 1 ). 

Option 1 provides an extension of the 
methodology, leading to the estimation 
of F, E (= F/Z) and probability of 
capture; here, an estimate of the mean 
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Functions 


annual habitat temperature, (T, in °C) 
will be required, to estimate M from 
Pauly's empirical formula (see further 
below). 

Option 1: Without seasonality 


ln(Nj/Atj) = a + b • tj 


where 

N is the number of fish in length class i, 
At is the time needed for the fish to 
grow through length class i, t is the age 
(or the relative age, computed with 
t 0 = 0) corresponding to the mid-length 
of class i, and where b, with sign 
changed, is an estimate of Z. 

Following estimation of Z, the routine 
can be used to estimate M using Pauly's 
M equation and F, from Z = M+F, as 
well as the exploitation ratio, E = F/Z. 

Catch curve analysis can then be 
extended to an estimation of 
probabilities of capture by backward 
projection of the number that would be 
expected if no selectivity had taken 
place (N'), using 


N'j_| = N'j • EXP(ZAt) 


with At as defined above and 


Z = (Zj + Z i+ j)/2, 

Zj = M + Fj, 

Fj_ i = Fj - X, and 

X = F / (no. of classes below P] +1); 
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Outputs 


User interface 


and where P| is the first length group 
with a probability of capture equal to 
1.0, and whose lower limit is an 
estimate of L'. From this, probabilities 
of capture by length are computed from 
the ratios of N/N'j. 

Option 2: With seasonality 


ln(N) = a + b • f 

where N is the number of fish in 
(pseudo) cohorts "sliced" by means of 
successive growth curves, t' the relative 
age of the fish in that pseudo-cohort, 
while b, with sign changed, provides an 
estimate of Z. 

Option 1 

Estimates of Z, F, M. E, L', and 
probabilities of capture. 

Option 2 

Estimate of Z 

The first tab of the length-converted 
catch curve routine (Fig. 4.20) is similar 
to that presented for Powell-Wetherall’s 
plot (Fig. 4.5) except that the growth 
parameters are required inputs. Option 
2 will be taken if parameters indicating 
seasonality in growth are given. 
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Fig. 4.20. The length-converted catch curve user interface in FiSAT II. 

The growth parameters are required to 
compute the data points presented in the 
second tab. 

If the required inputs are encoded, 
clicking the “Compute Data Points” 
will plot the results presented in the 
second tab (Fig. 4.21). 
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I V Length-converted Catch Curve [C: ' Program 




Length-Converted Catch Curve 


Header 

Filename: |CAProgram FiesVFtSAT ll\Data\Sampte Iq 

Specie* name: jUpeneus moiuccemjs 

Other identifiers: |RAGAY GULF 1980 


*1 



Fig. 4.21. Length-converted catch curve based on the growth parameters entered in the 
first tab of the user interface. 


The first and last data points to be 
included in the regression analysis are 
computed automatically but users have 
the option to change the default points 
by clicking on the “Reset Selections” 
command button. The first and last data 
points may be selected or deselected 
using the mouse. 

The statistical results of the regression 
analysis are presented in the third tab. 
When C is not given or set to zero, the 
regression line can be extrapolated to 
approximate the probability of capture 
given natural mortality. FiSAT II 
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provides the option to estimate this 
value using Pauly’s empirical equation 
(Fig. 4.22). 



Fig. 4.22. User interface to compute natural mortality (M) 
from Pauly's empirical equation. Users may encode the 
value of M directly that will be used in extrapolating points of 
the catch curve (Fig. 4.23) 

The extrapolated points (Fig. 4.23) will 
be used to approximate the probability 
of capture. To estimate the selection 
parameters based on the extrapolated 
points, click the command button 
“Show Prob.” that replaced the 
command button “Extrapolate Prob.” 

The probability of capture for each 
length group will be plotted and users 
have options to use either the logits plot 
or simply using the running average 
technique to estimate selection 
parameters (see below for detailed 
discussion on how to use the user 
interface in estimating the selection 
parameters). 
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I rngth -converted Catch Curve [C:\Program F 


JSSSESSZ 

Length- Converted Catch Curve 


Xj 


Header 



Filename. 

jCAProgram Fies\FiSAT ll\Data\Sample Jq 


Species name: 

|Upeneus mduccensis 


Of her identifiers 

|RAGAYGULF 1980 


Qata and Options 

Qjtch Curve j Regression Statistics | 



Length-Converted Catch Curve 
(lor Z=6 56, M (at 27 0*0-1 .94. F=4 62, E=0 70) 


Growth Parameters 
Loo : 



(L 

C: 

WP: 

to: 


Pont Parameters 


1st 


□ass no 
Y: ln(N/Zit) 

X relative age 


880 


Last 


13 

3.73 

179 


1 03 

Reset Selections 


Z Estimate 

Z I 656 

□ olZ | 480-| slT 


mpute Oata Poi 


£howProb 


Erint.. 


Close 


Fig. 4.23. Length-converted catch curve with extrapolated data points. Clicking "Show 
Prob." will open the dialog box to estimate the selection parameters (see further below 
for details). 


Remarks In both options, defaults for the points 

to be included in the regression analysis 
are proposed by FiSAT II; use these 
only if no better sequence of points can 
be identified by visual inspection of the 
catch curve. 


Option 1 provides the facility to 
proceed to the estimation of L 2 5 , L 50 
and L 75 . This routine is discussed 
further below. 

Readings Pauly (1984a; 1984b; 1990) 




95 


Jones/van Zalinge plot 

The cumulative plot of Jones and van Zalinge is an early 
form of length-converted catch curve and these two 
methods thus share many common assumptions. 

Required file Length-frequency data with constant 

class interval. 


Input parameters L, and K (or 1 if K is unknown), with 
graphical selection of the points to be 
used in the regression analysis. 


Function 


ln(C Li 00 ) = a + b • In (L.-Lj) 


where 


C^j x is the cumulative catch (computed 
from the highest length class with non- 
zero catch) corresponding to length 
class i, 

and 

Lj is the lower limit of length class i. 
The slope b, is an estimate of Z/K. 

Outputs Z or Z/K if 1 was entered instead of K. 

User interface The first tab of this user interface 
identifies the file and provides options 
(see Fig. 4.5) for temporary data 
transformation. 

The data points will be plotted after 
encoding the growth parameters to be 
used in the analysis. Points to be 
included in the regression analysis are 
not automatically identified by FiSAT 
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II as in the case of the length-converted 
catch curve routine. The user will be 
required to identify the first and last 
points to be included using the mouse 
(Fig. 4.24). 

The details of the regression analysis 
are given in the third tab. 



Fig. 4.24. User interface of the Jones and van Zalinge plot to estimate total 
mortality. 

Remarks This routine was included in FiSAT II 

mainly for comparison, as it lacks many 
of the advantages of length-converted 
catch curves (independence of the data 
points, ability to deal with seasonal 
growth oscillations, etc.), without 
offering any distinct advantage over 
catch curves. 


Readings Jones (1984) 

Jones and van Zalinge (1981) 
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Z from mean length 

Two routines are available in FiSAT II which derive 
estimates of Z from mean lengths: 

Beverton and Holt model 

This model, a classic developed by Beverton and Holt 
(1956), assumes that growth follows the VBGF, that 
mortality can be represented by negative exponential decay, 
and that L is estimated from a sample representing a 
steady-state population. 

Required file None 

Input parameters Growth parameters L x and K, cut-off 
length (L’) and mean length (L) above 
L'. 

Function Z = K (L,-L)/(L,,-L') 

Output Estimate of Z 

User interface This routine will require only the inputs 

to the growth parameters (L and K), 
mean length and the cut-off length to 
estimate the total mortality. Note that if 
K is unknown and set to 1, the result 
should be read as Z/K. 

If a value is altered, clicking the 
“Compute” command button will re- 
compute the Z estimate (Fig. 4.25). 
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Fig. 4.25. User interface in estimating the total mortality 
using the Beverton and Holt model. 

Remarks This routine is presented along with the 

model of Ault and Ehrhardt (see below) 
because the inputs are similar. This 
model is best applied to long living and 
slow growing species. 

Reading Beverton and Holt ( 1 956) 

Ault and Ehrhardt method 


Unlike the Beverton and Holt Z equation (above), this 
model does not assume an infinite life span for the fish of 
the stock being analysed and thus is applicable to short- 
lived tropical species. 

Required file None 

Input parameters Growth parameters L x and K, the cut- 
off length (L’X the mean length (L) and 
the maximum length (L max ) 
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Functions 



where 

A(L') = Z(L'-L)+K(U-L), and 
A(L max ) = Z(L max -L)+K(L x -L) 

Outputs Estimate of Z 

User interface As with the previous routines, the total 
mortality will be estimated only after 
the required inputs are given (Fig.4.26). 


■gff.M-WW xj 


Related to the B evert on and Ho* model, n the Ault and 
Ehrhardt model. Z from mean length, is 


Loo l««i 

Z / K 

Zff L*e«n| ♦ KfLoo Lmean) 

Loo - L* 

Z[[ man Lmean) ♦ K.|Loo l aiean 


where Loo «s the asymptote length. K is the cur vatu e 
parameter of the VBGF, Lmean « the mean length ol the fish m 
a sample representing a steady-state population. L' it the cut-off 
length or the lower fcmit of the smalest length clast included in 
the computation, and Lmax is the largest fish in the sample 


Parameters 



Asymptotic length (Loo) 

1 

i 

Growth constant (K/year) 

1 


Mean length (Lmean) 

1 


Cut-off length (Lf 

r~ 


Maximum length (Lmax) 

i 


Computed total mortatty £) [" 


Cgmpute | Close 


Fig. 4.26. User interface to estimate the total mortality 
using the Ault and Ehrhardt model in FiSAT II. 
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Remarks The function is solved iteratively. The 

value of L max required here may be the 
largest specimen in a set of samples, or 
L’max as estimated via the extreme 
value theorem (see Chapter 5). The 
required L and L' can be obtained as for 
the Beverton and Holt model (see 
above). 

Readings Ault and Ehrhardt (1991) 

Ehrhardt and Ault ( 1 992) 

Hoenig's model 

The method below relates t max (i.e. the longevity of fish) 
and mortality. 

Required fde None 

Input parameters t max (maximum age in years) observed 
in a given stock and as an option, t c 
(mean age at first capture) and the 
sample size, N from which t max was 
estimated. 

Functions When t c and N are not available, Z is 

computed from the empirical model 

ln(Z)= 1.44 - 0.984 ln^) 

When t c and N are available, total 
mortality is estimated from 

z=i/(ci-(t max -g) 

and the standard error (s.e.) of the 
estimated Z is from the relationship: 


s.e.(Z) = (c 2 -Z 2 )'/’ 
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where ci and C 2 are table values and a 
function of N (these table values are 
pre-programmed in FiSAT II). 

Output Estimate of Z 

User interface The user interface of this routine 
provides the user with options on what 
model to use as presented by the author 
of the model. Note that the text boxes 
for N and tc will be enabled only if the 
checkbox for Model 2 is checked using 
the mouse (Fig. 4.27). 


Remarks The estimates of Z obtained by these 

models (particularly the first) are rather 
approximate, and should be validated 
using other methods. 



Fig. 4.27. User interface of Hoenig's model as 
implemented in FiSAT li. Note that the estimates using 
the two models can be calculated simultaneously. 
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Tabulated values of C| and C 2 may be 
found in Hoenig and Lawing ( 1 982) and 
Pauly (1984a). 

Readings Hoenig ( 1 982) 

Hoenig and Lawing ( 1 982) 

Pauly (1984a) 


Natural mortality 

In fish stock assessment, natural mortality is the most 
difficult parameter to estimate. The following models, two 
empirical and one analytical, may be useful to 
approximate M. 

Rikhter and Efanov’s method 


A variant of the Beverton and Holt (1956) equation which 
relates longevity (t max ) with the L mass /L x ratio (where 
Lmass ' s the mean length at first maturity), this model 
relates natural mortality (M) to the age at which 50% of the 
stock reaches the age of "massive spawning" (t mass ). 

Required file None 


Input parameter t mass (*= age at first maturity) 

Function M = ((1.52/t mass ) 0.72) - 0. 16 

Output Estimate of M 


User interface The dialog box for this routine will 
require only the input of the age at 
massive maturation to estimate the 
natural mortality (Fig. 4.28). 
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Fig. 4.28. User interface to estimate natural mortality using 
Rikhter and Efanov's empirical model. 

Remarks The derivation of the model was based 

largely on data from boreal/temperate 
stocks and should not be used for 
tropical stocks without prior 
comparison with the results of other 
models. 

Reading Rikhter and Efanov ( 1 976) 

Pauly's M equation 

These equations were derived from 175 independent sets of 
estimates of M and predictor variables for most tropical 
species. Two options are provided; (1) for cases when L. is 
available and (2) for cases when W* is available. 

Required fde None 

Input parameters Growth parameters L„ (in cm) or W, 
(in g) and K (year 1 ) and T [mean 
annual habitat temperature (in °C)]. The 
values of L. to be used should refer to 
or approximate total length. 
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Functions Option 1: L* is given 

ln(M) = -0.0152 - 0.279 ln(L x ) + 

0.6543 ln(K) + 0.463 ln(T), 

and 

Option 2: W x is given 

ln(M) = -0.485 1 - 0.0824 ln(W x ) + 

0.6757 ln(K) + 0.4687 ln(T) 

Output Estimate of M 

User interface The user interface of this routine 
provides the users with the two options 
to estimate natural mortality (M) based 
on Pauly’s M empirical equations. 
Option 1 is the default model. Clicking 
on the second option (red circle in 
Fig. 4.29) will change the labels of the 
function and input fields. 



Fig. 4.29. User interface to estimate the natural mortality using Pauly's M 
empirical equations. 
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Remarks These equations provide reasonable 

estimates of M for fish, and perhaps 
shrimps and squids as well, but should 
not be applied to bivalves for which 
M»K may provide a better 
approximation. The models assume L* 
to be in cm, and W« in g; FiSAT II 
automatically performs any required 
transformation. Also, values of T down 
to -2°C can be entered and are adjusted 
internally. 

Readings Pauly ( 1 980, 1 984a) 

Recruitment patterns 

This routine reconstructs the recruitment pulses from a time 

series of length-frequency data to determine the number of 

pulses per year and the relative strength of each pulse. 

Required file Length-frequency data with constant 

class size. 

Input parameters Growth parameters: L,, K, (C, WP), and 
t 0 if available. 

Functions Backward projection, along a trajectory 

defined by the VBGF, of the 
frequencies onto the time axis of a time- 
series of samples (one single, but well- 
structured sample may also be used). 

Two options are provided; 

Option 1 , use of the length-frequency 
data. 

Option 2, use of restructured data. 
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Output Plot showing the seasonal pattern of 

recruitment. 

User interface The user interface for this routine 
contains three tabs. The first tab is to 
identify the file to analyse and options 
for data transformation (Fig. 4.30). 
Another option is provided to the user, 
i.e. the use of the length frequency data 
as it is and use of restructured data as is 
done in ELEFAN I. 


? » Recruitment Pattern... 


Recruitment Pattern 

£ata and Option: j Growth Parameter Inputs & Graph | Numeric RenJt: j 


*J 


Header 



Ftename 

|c AProgram Fies\FiSAT ll\Data\Sample Kq 


Species name 

|Upeneus mohjccensts 


Othei identifier: 

|RAGAY GULF 1980 



Sample Selection 

S electy deselect data to ndude m the 
analysts 


</ 2/3/1980 
✓ 21/5/1980 
«> 23/7/1980 
f 27/9/1980 
•/ 20/11/1980 


Data Source 
(• Use original data 
/'Use restructured data 




Data Transformation 

(• Raw data (no transformation) 

(Use this when samples represent catch per 
effort) 

C Percent ot sample total 

(Use this when samples have very dissimiar 
sizes) 

C Square-root of sample total 

(Use this when samples have roughly srmiar 
sizes) 

f~ Apply temporal weiglhng ?| 


Compute | Print. | Close 


Fig. 4.30. File identification tab and options to temporarily transform data before data 
analysis. 


The growth parameters are required 
inputs and clicking the “Compute” 
command button will compute the 
recruitment pulse(s) (Fig. 4.31). In 
cases where there are two recruitment 


Copyrighted malirial 




107 


pulses, the composite data can be 
decomposed using the approach as 
implemented in NORMSEP (see 
above). In this case, starting estimates 
of the mean (in months) are required to 
compute the Gaussian distributions. 
The numeric results are given on the 
third tab. 

When t 0 is set equal to zero, the 
ordinate scale of the plot will be 
relative, not calendar time. 



Fig. 4.31. User interface to estimate the recruitment pulses in FiSAT II. 


Remarks The results generated by this routine 

should be treated as approximations 
because, while allowing statements on 
the number of annual pulses and on 
their relative strength, this model is 
based on two assumptions that are 
rarely met in reality: (i) all fish in the 
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sample grow as described by a single 
set of growth parameters and (ii) one 
month out of twelve always has zero 
recruitment. When using restructured 
data, the temporal spread reduces and 
thus probably better reflects the actual 
seasonality of recruitment. As an 
additional option, FiSAT II also allows 
the user to analyse recruitment patterns 
using a maximum likelihood approach 
(see section on NORMSEP) to fit a 
Gaussian distribution. 

Readings Moreau and C uende (1991) 

Pauly (1983) 

Probabilities of capture 

"Probabilities of capture..." leads, in FiSAT II, to those 
routines that can be used to determine gear-specific 
selection curves. 

T rawl-type selection 

Trawl-type selection curves are characterized by left-hand 
selection, i.e. small fishes escape through the mesh and 
hence are not caught. In most cases, this type of data is 
obtained from trawl selection experiments, which usually 
consist of covering the cod-end whose selectivity is to be 
assessed with a fine-mesh cover. However, similar 
information may also be obtained from length-converted 
catch curves, by backward extrapolation of the catch curve, 
and comparison of the numbers actually caught with those 
that "ought" to have been caught. 

To estimate selection parameters, (L 25 , L 50 , L 75 ), FiSAT II 
provides the user with the options of either: ( 1 ) using the 
logistic curve, which assumes selection to be symmetrical 
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or nearly so or (2) using moving averages over three classes 
and interpolating the selection parameters. 

Required file Probability data file 

Input parameters Option 1: Logistic curve 

Visual selection of the points of a logits 
plot to be used for the regression 
analysis. 

Option 2: Moving average 
None 

Option 1: Logistic curve 

ln((l/P L )-l) = Sl -S2-L 

where P L is the probability of capture 
for length L, and 

L 25 = (ln(3)-Sl)/S2 
L 50 = S1/S2 
L 75 = (ln(3)+Sl)/S2 

Option 2: Moving average 

^L,i(new) — (^L.i-l^^L.itoldl^^L.i+l)^’ 

which provides a smoothed series of 
probabilities from which L 2 5 ? L 5 q and 
L 75 are estimated through linear 
interpolation. 

Estimates of L25, L50 and L 75 

This routine (Fig. 4.32) requires a 
probability of capture file that can be 
created based on selection experiments 


Functions 


Outputs 
User interface 
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or from a length-converted catch curve 
routine (see above). This routine can 
also be accessed directly from a length- 
converted catch curve routine. 

Two options are provided to the user to 
estimate the selection parameters. The 
use of the logistic transformation plot 
(also referred to as logits; Fig. 4.33) 
will require the user to identify at least 
three data points to be included in the 
analysis (default is all data points will 
be used). The use of the running 
average technique will not require 
additional user intervention to estimate 
the parameters. 


PtotMMrtp ot L<sp(uic 


*J 


H«ade» 

tfc n a w |r ll\Dtfa\S«npie prt> 

Spec*! name jupeneu* mokxxmm 



I Pw I 


Fig. 4.32. The user interface for the estimation of the selection 
parameters using either logistic transformation or running 
average technique. 
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Fig. 4.33. Logistic transformation of the data points to 
identify points to be included in the analysis. 

Remarks This routine may also be accessed 

directly from the length-converted catch 
curve analysis. The selection curve can 
be saved to a file for subsequent 
analysis of L/F data for the effects of 
gear selection and incomplete 
recruitment or as an input to yield per 
recruit analysis. 

Reading Pauly ( 1 984a) 

Gillnet selection 

The selection curves of gillnets (and some other gears such 
as hooks) are characterized by strong left-hand and right- 
hand selection. This selection type can be modelled using 
normal (Option 1) or log-normal (Option 2) selection 
curves. 

Required file Length-frequency data file with 

constant class size with two samples 
representing the catches (in numbers) of 
two gillnets of different mesh sizes, 
with at least 3 overlapping non-zero 
length classes. 
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Input parameters 


Functions 


Mesh sizes of the two selected samples 
(m A and m B ), where one (m A ) should 
be smaller than the other. 

The probability of capture for mesh size 
m A and m B at a given length L is, 

P A, L = ex P (-((L-L A ) 2 /(2a 2 )), 

P B, L = ex P (-((L-L B ) 2 /(2 a 2 )), 

where L A and L B (the optimum length 
for mesh sizes m A and m B respectively) 
are given by 


L A = SF • m A , 
L B = SFm B , 


and where SF is the selection factor 
computed from, 

SF = (-2a)/(b(m A +m B )). 

The coefficients a and b are estimated 
from the regression: 


ln(Cj B /C i A ) = a + b • Lj 

where the index i denotes length 
classes, and Cj A and C; B are the 
observed catches (in numbers) for class 
i of gears A and B, respectively. 

The standard deviation (a), assumed 
equal for both mesh sizes, is then 
computed from 

o = ((2a(m A -m B )/(b(m A +m B ))) 1 ' 2 
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For Option 2, the logarithms of the 
variables referring to length (Lj, m A , 
mg) are first taken before the functions 
are applied. 

Outputs Probabilities of capture by length for 

each mesh size and the joint standard 
deviation (a). 

User interface The user interface of this routine has 
two tabs. The first tab (Fig. 4.34) is to 
identify the length frequency file (with 
at least two samples). In addition, the 
two samples have to be identified and 
corresponding mesh sizes encoded. 


JSi 


Gillnet Selection 



Cflwouig | girt | 5*»e | Cbw 


Fig. 4.34. User interface to identify the length frequency file, the samples and 
corresponding mesh sizes used to generate the catches. 
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Once the required inputs have been 
entered, clicking on the “Compute” 
command button will plot the 

distribution (Fig. 4.35) and estimate the 
probability curve that can be saved as 
probability data for further analysis. 



Fig. 4.35. Computed probability curve based on catch data from two gillnets of 
different sizes. 

Remarks The mesh sizes (m A , m B ) will be read 

from the "Gear description" field in 
sample sub-headers (default values) if 
they are not entered at run time. Also, 
note that the normal and log-normal 
models are performed concurrently and 
that the user can toggle between the 
results from both models. 


Readings Pauly ( 1 984a) 

Pauly (1983) 


I 
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Virtual population analysis 

Virtual population analyses (VPA) are methods which 

allow the reconstruction of the population from total catch 

data by age or size. 

Age-structured VPA 

This is the version of VPA proposed by Gulland (1965). 

Required file Catch-at-age data file 

Input parameters Natural mortality (M), an initial guess 
of terminal fishing mortality (FJ, and 
the time interval between age groups (in 
years). 

Functions The reconstruction of the population 

starts by estimating the terminal 
population (N t ) given the inputs, from 


N t = C t • (M + Ft)/F t 


where C t is the terminal catch (i.e. the 
catch taken from the oldest age group). 

Then, starting from N t , successive 
values of F in the previous age group 
are estimated, by iteratively solving 

q = N i+1 • (Fj/Zj) • (EXP(Zj)-l) 

with population sizes (Nj) computed 
from. 


Nj = N i+1 • EXP(Zj) 


The last two equations are used 
alternatively, until the population sizes 
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and fishing mortality for all age groups 
have been computed. 

Outputs Plot of the reconstructed population and 

fishing mortality for each age group. 

User interface The routine (Fig. 4.36) requires catch- 
at-age data and inputs of the required 
parameters (in years). The second tab 
(Fig. 4.37) presents the results of the 
analysis. The stacked histograms 
present the survivors, population lost to 
natural causes and the catches. Overlaid 
in the stacked histograms is the 
computed F-array. The numeric results 
are presented in the third tab. 


Age-structured VPA. 


Age-Structured Virtual Population Analysis (VPA) 

Data | Graphic Results | Tabular Results | 




neaaer 



Ftename 

|CAProgram Fdes\FiSAT ll\Oata\S ample caa 


Species name 

|Sardna pdchardus 


Other identifiers 

| Caught off Morocco 


Umt ol catch measuement 

|nurrfc»s Und ot age meijuiement [ ye «i 



Parameter 


Natural mortality (M /ytl | 8 

T eimmal fishng mortality (Ft /yv): | .3 

T roe interval between age groups (n yr) P 


.25 


Proceed 


Pnnt.. 


£lose 


Fig. 4.36. File identification and input of parameters when analysing catch-at-age 
data using the age-structure VPA in FiSAT II. 
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Fig. 4.37. VPA plot of the results. Note the option to display the legend and 
position it anywhere on the graph. 

Readings Gulland (1965) 

Mesnil (1985) 

Pauly (1984a) 

Pope (1972) 

Length-structured VPA 

This routine modified from Jones and van Zalinge (1981) 
utilizes basically the same approach as the previous routine 
(age-structured VPA), but is adapted to accommodate 
length frequencies. 

Required file Length- frequency data file (representing 

mean annual catch at length, see 
below). 

Input parameters L, and K. Note that L, must be at least 
10% larger than the largest fish in the 
file. If you need to increase L„, use o' to 
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Functions 


Outputs 


reduce K. accordingly. Alternatively a 
plus group may be created. 

The initial step is to estimate the 
terminal population (N t ) given the 
inputs, from 

N t = C t • (M + F t )/F ( 

where C t is the terminal catch (i.e. the 
catch taken from the largest length 
class). 

Then, starting from N t , successive 
values of F are estimated, by iteratively 
solving, 

Cj = N i+At • (Fj/Zj) • (EXP (Z; Atj)-1) 
where 


Atj = (t i+1 - tj), 


and 

tj — t 0 - (1/K) ■ ln( 1 -(L/L*)) 

where population sizes (Nj) are 
computed from 

N i = N i+At EXP(Z i ) 

The last two equations are used 
alternatively, until the population sizes 
and fishing mortality for all length 
groups have been computed. 

An F-array representing the fishing 
mortality for each length group, the 
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reconstructed population (in numbers), 
and the mean stock biomass by length 
class. 

User interface The user interface of this routine 
contains three tabs. The first tab (Fig. 
4.38) is to identify the file to use and 
data transformations. Note that if data 
were raised by some value, this should 
be re-transformed to represent catches 
(in numbers; red circle in Fig. 4.38). 


*1 


Length-Structured Virtual Population Analysis (VPA) 


Header 
Filename 
Species name 


|c \Program FilesVFiSAT ll\Data\Sample Hq 
|Upeneus moluccensis 


Other identifier*: |RAGAY GULF 1980 

£>«ta and Options | Parameters l VPA Plot] Jabular Results | 


Data Selection 

WARNING! The catches |m numbers) 
should represent the mean total annual 
catch 

The catches (n numbcis) are tn 
1 * 10 A | 6 

Select or deselect samples to include in 
the analysis 

* 20/1/1980 
v» 2/3/1980 

✓ 21/5/1980 
23/7/1980 

v 27/9/1980 

✓ 20/11/1980 

vi 


Data T rans/ormation 

6 Raw data (no transformation) 

(Use this when samples represent catch per 
effort) 

r Percent of sample total 

(Use this v4ien samples have very dtttimiaf 
sees) 

f Square-root of sample total 

(Use this when samples have roughly svniiar 
sees) 

r Apply temporal weigthng ?| 


Cflmpute 


| p°« 


Fig. 4.38. Data identification tab of the length-structured VPA routine in FiSAT II. 


The required parameters, including the initial 
estimate of the terminal fishing mortality 
(Ft) are given in the second tab (Fig. 4.39). 
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As in the age-structured VP A, the legend can 
be displayed and re-positioned on the graph. 


To alter the value of the terminal F 
graphically, check the “Modify terminal F 
from Graph” option and click on the graph to 
represent the position of the terminal F. If 
the mouse is pointing to a length group less 
than the largest recorded length, a +Group 
will be generated. 


*J 


Length-Structured Virtual Population Analysis (VPA) 


|c AProgram FilesSFiSAT ll\Data\Sample Iq 


Header 
Filename 

Species name |Upeneus motaccenas 
Other »dentrf»ers IRAGAY GULF 1980 


P Display legend V Modrfy terminal F from graph 





Enter the 
coefficients of the 
length-weight 
relationship [W ■ 
aL A b* 



Cgmpute J $ave Results | Print 


Close 


Fig. 4.39. Results of the length-structured Virtual Population Analysis (VPA). 
The terminal F may be modified directly from the graph. 


The numeric equivalent of the graph is 
given in the third tab. 
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The resulting F-array can be saved as a 
Species Table that can be accessed 
directly when creating Scenario Files 
for the Thompson and Bell routine of 
FiSAT II. 

Remarks The "length-frequency" file used here 

must in fact consist of total annual 
catch-at-length data, in numbers. It may 
be the average total catch of several 
years. 

Reading Pauly ( 1 984a) 


Predictions 

While the previous routines of ASSESS are used to 
estimate the value of certain parameters, the routines 
below, which require these inputs, are meant to be used for 
yield and stock predictions, and hence, to identify 
appropriate management regimes. 

Relative Y/R and B/R analysis: Knife-edge selection 

The relative yield-per-recruit model based on the Beverton 
and Holt model of 1966. 

Required file None 

Input parameters L c /L, ratio (from 0.05 to 0.95) and M/K 
ratio (from 0.10 to 9.99) 
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Functions 


Y'/R = 


Outputs 


User interface 


Relative yield-per-recruit (Y'/R) is 
computed from: 

(l + m) (l + 2m) (l + 3m)| 


where 


U= 1-(L C /L.) 
m = ( 1-E)/(M/K) = (K/Z) 

E=F/Z 

Relative biomass-per-recruit (B'/R) is 
estimated from the relationship 

B'/R = (Y'/R)/F, while 

E max , E 0 ] and E 0 5 are estimated by 
using the first derivative of this 
function. 


Plots of Y'/R vs E (=F/Z) and of B'/R vs 
E, from which E max (exploitation rate 
which produces maximum yield), E 0 | 
(exploitation rate at which the marginal 
increase of relative yield-per-recruit is 
1/1 0 th of its value at E=0) and E 0 5 
(value of E under which the stock has 
been reduced to 50% of its unexploited 
biomass) are also estimated. 

FiSAT II automatically plots the Y’/R 
isopleths diagram with M/K set to 1 .00. 
If the user changes the M/K value (red 
circle in Fig. 4.40), a new isopleths 
diagram will be plotted automatically. 
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Fig. 4.40. Y7R isopleths diagram with M/K value set to 1.0 (see figures 4.41 and 
4,44 for facsimile presentations of other possible options In the presentation of 
the results). 


The user interface also provides other 
options for viewing the results. Fig. 4.41 
is an example of a 3D presentation of 
results for B’/R using the same M/K 
value. Note that the ceiling display, 
floor display and spin plot options are 
applicable only when the plot is 
presented in 3D. 

To manually spin the plot in all axes, 
click both buttons of the mouse and drag 
the plot to the desired viewing angle. 

Whether the results are presented as 
isopleths or in 3D, moving the mouse 
over the plot identifies the numeric 
values of the plot. 
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A 2D plot of the results (see a similar 
plot in Fig. 4.44) may also be plotted by 
clicking on the “2D Analysis” command 
button. 



Fig. 4.41. 3D presentation of the B'/R plot with a floor display of the results. 

Remarks An example of the 2D graph is given in 

Fig. 4.44. 

Readings Beverton and Holt (1966) 
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Relative Y/R and B/R analysis: Using selection ogive 

The relative yield-per-recruit model presented in the 
following is based on the Beverton and Holt model of 
1966, modified by Pauly and Soriano (1986). 

Required file Probabilities of capture data 

Input parameters L«, K and M 

Functions Relative yield-per-recruit (Y'/R) is 

computed from 


Y'/R = 2P i ((Y7R) i -G i _ 1 M(Y'/R) i+1 -G i )) 


where 


(Y’/R)j refers to the relative yield-per- 
recruit computed from the lower limit 
of class i using 


(Y'/R), = EU M/K 


3U 3U 2 U 3 j 
(l + m ) (l + 2m) (l + 3m)j 


where U and m are defined as above, 

Pj is the probability of capture between 
Lj and Lj + i, while Gj is defined by 

G; = n rj 

where 

rj = ( 1 -Cj) Si /( 1 -C|_ | ) Si , and 
Sj = (M/K)(E/(l-E))Pj. 

Here, B'/R is estimated from 
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(B'/R)j = (l-E)A/B 
where 


A = 


3U | 3U 2 U 3 1 
(l + m) (l + 2m) (l + 3m)| 


R_fi 3U 3U : U J 1 

| (l + m') (l + 2m') (l + 3m')| 

and where m' = 1/(M/K) = m/(l-E). 

E max , E 0 i and E 05 are estimated by 
using the first derivative of the function. 

Outputs Plots of Y'/R vs E (=F/Z) and of B'/R vs 

E, from which the E max , E 0 [ and E 0 5 
(as defined above) are also estimated. 

In addition to these two-dimensional 
outputs, the user has the option to plot 
yield-isopleths diagrams, which can be 
used to assess the impacts on yields of 
changes of E (corresponding to the level 
of exploitation) and of c (- L 50 /L., 
corresponding to a change of mesh 
size). 

User interface The user interface of this module is very 
similar to the presentation of the results 
using the knife-edge assumption. The 
display options are the same. However, 
a selection data (probability tile) has to 
be loaded before any calculations are 
done by the routine (Fig. 4.42). 
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Fig. 4.42. Initial display of FiSAT II when the Beverton and Holt Y/R model using 
selection ogive routine is accessed. A probability of capture file is required by the 
routine. 


The L. used in previous calculations 
using the file or associated files are used 
as defaults. The user may change the 
value to compute new points of the plot. 

As in the previous routine, the results 
may be plotted in three dimensions (Fig. 
4.43) with options to also plot the 
ceiling and floor displays. 

The 2D plot (Fig. 4.44) is the graphical 
representation of the results for a 
specific c (= L 5 o/L*) value. The user 
has the option to change the value or 
may slide the scroll bar to examine 
changes of the c value (simulating a 
change of the mesh size). 
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Fig. 4.43. A 3D presentation of the Y'/R plot using the selection data. 



Fig. 4.44. 2D presentation of the results for specific Lm/L. value. The 
scroll bar may be used to examine changes of the Lm/L. values. 
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Remarks Using probabilities of capture which 

abruptly change from zero to 1 at 
L 50 (=L c ) enables this model to simulate 
the behaviour of the same model 
assuming a knife-edge selection, and 
thus allows comparison between the 
two approaches. 

Readings Pauly ( 1 984a) 

Pauly and Soriano (1986) 

Beverton and Holt (1966) 

Silvestre et al. (1991) 


Thompson and Bell yield and stock prediction 

This model combines features of Beverton and Holt's Y'/R 
model with those of VPA, which it inverts. The version 
presented here can be used to analyse either a single 
species, exploited by a single gear, or several species, 
exploited by several fleets. Naturally, the data requirements 
increase with the complexity of the analysis required. It is 
for this reason that this routine is presented through a series 
of "options". 

Five options are provided: (1) Create, edit, save or print a 
Species Table, containing population parameters, value by 
length groups and the F-array associated with one species, 
(2) Create, edit, save or print a Fleet Table, containing a 
text description of the fleet, (3) Create, edit, save or print a 
Relations Table linking the Species and Fleet Tables, data 
on selection and catch indexes used, to split the F-array if 
the species is (are) exploited by several fleets, (4) Run a 
Predict routine, which executes the model and outputs the 
graphical and numeric results, (5) Help, which contains a 
short description of the procedure and data requirements. 

Required file Scenario File 
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Input parameters 


Functions 


Species Table(s), Fleet Table(s), and 
relations between the species and the 
fleets. 

The sum of the yields (Y = 2Yj) is 
computed from 

Yj = Cj • Wi 

where the mean body weight 



and where a and b are the coefficients 
of the length-weight relationship and Lj 
and L j+ ] are the lower limit and upper 
limit of the length class, respectively; 
also we have 

C i = (N i -N i+1 )(F i /(M + F i )) 

where the predicted population (N ; ) is 
given by 


Nj +] = Nj • EXP(-(M+Fj)-Atj), and 


Atj = (1/K) • ln((L x -Lj)/(L*-Lj +1 )) 


The biomass is computed from 

B i =((N i -N i+1 )/(M+F i ))-At i w I 

and 

the value (V;) is computed by 


Copyrighted material 



131 


Outputs 


User interface 


V = Y- • v- 

v i 1 1 v i 

where Vj is the unit value for class i. 

In the multi-species/multi-fleet 
scenario, an F-array has to be generated 
(based on the ratio of catches of the 
different gears) before any yield 
predictions can be made. The yields, 
computed on a per species and per fleet 
basis, are then added. 

Plots of yields, values and biomass 
estimates for a range of f-factors for 
each species-fleet combination as well 
as the cumulative curve. An option is 
provided to examine the effect of 
changes of selection parameters. 

The user interface of this routine has 
three tabs. The first tab (Fig. 4.45) 
opens the scenario file (*.SCE) that 
identifies the list of species and fleets 
included for analysis. 
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llionipuMi imnI Hrll Y-S PrcAcIton.- 


*J 


Thompson and Bell Yield - Stock Prediction 
General | Gtaphcai Rents] Iabtiar Rents) 


bcenano 



Fiename 

C \Piogram FlesNFiSAT ll\Oata\$arrple tee 


Tide 

M iti-Speae* fleet 


Remarks 

This is orly a umJafed He 

3 



3 


Spaces Lstl 


Red Ld 


Iiarch T rad Fled (vay Mac»or| 


Compute ^rnitfe Sdectrnty Pirl Fie Qeabon Help Pose 


Fig. 4.45. User interface of the Thompson and Bell yield-stock prediction module that 
identifies the file and options to simulate changes in the f-factor for each fleet included in 
the scenario. 


As an option, the user may also uncheck 
the box corresponding to the fleet label 
to simulate a fixed value for the f- 
factor. In which case, the user will be 
prompted on what value to use. If the 
box is checked, FiSAT II will present 
results with the f-factor of the fleet 
varying from 0 to 4. 

When the “Compute” command button 
is clicked, the results will be plotted as 
shown in Fig. 4.46. The numeric 
equivalent of the result is given on the 
third tab. 
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Fig. 4.46. The graphical result of the Thompson and Bell yield-stock prediction module of 
FiSAT II. The left panel with three graphs shows results specific to a fleet and a species 
while the right panel shows the cumulative results. 


The three graphs in the left panel are 
the results of the simulation for a 
specific species as exploited by a fleet. 
Use the dropdown list in the upper left 
panel of the interface to view results of 
other combinations. The range of 
colours used (red to light blue; the blue 
band is not visible in Fig. 4.46) 
represents the age group distribution 
with red representing the oldest age 
group. 

The right panel is the graphical 
representation of the cumulative results. 
The important points of the plotted lines 
are given by small circles. The small 
circle in the value and yield curves 


Copyrighted material 


Remarks 


Readings 


i 


134 


represents the maximum point while 
that for the biomass curve, represents 
the point at which the biomass has 
reach 50% of the original biomass. 

When used here, the output of a length- 
structured VPA fulfils most of the data 
requirements (on a per-species basis). 
Thus, this form of VPA was 
programmed to automatically generate 
the appropriate file for transfer to this 
routine. 

Sparre and Wiilmann (1993) 

Thompson and Bell (1934) 
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Chapter 5. SUPPORT menu 


What you will learn from this chapter 

This chapter will introduce you to routines meant to 
supplement or support the data analyses presented in the 
previous chapter. 

Simulation of length-frequency samples 

This routine applies the Monte Carlo technique to simulate 
the dynamics of a fish stock and a random sampling 
procedure, taking into account the various stochastic 
elements of the biological system when generating length- 
frequency samples. The file(s) output by this simulation 
routine may be used in various ways such as for assessing 
the applicability of a model to data with known 
characteristics or estimating population parameters by 
changing them iteratively until the length frequencies that 
are generated come very close to those observed. 

This routine includes options to generate length frequencies 
through random sampling in a "closed habitat" and length 
frequencies accounting for migrations. 

Required file None 

Input parameters Except for one element, the migration 
route for migratory stocks, the required 
inputs are the same for both options. 
These inputs are summarized in 
Tables 5.1 to 5.3. 

Functions The sequence of probability functions 

used are summarized in Table 5.4. 
These are used for simulating the life 
(growth, migration) of one individual 
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fish at a time, from one sampling period 
to the next, until it dies due to M or F. 

Tabic 5.1 . Required single field inputs for simulating length- frequency data. 


Field name 

Type 

Limits 

Remarks 

Number of areas 

N 

1 to 5 

The number of hypothetical areas 
between which the stock will 
migrate (for migratory stocks only). 

Sample type 

Ch 


Three sample types are available; 
(1) representing population, (2) 
representing catch, and (3) double 
sampling, representing both 

population and catch. 

Number of fish sampled 

N 

1 tolO ft 

Total number of specimens in a 
simulated sample. 

Number of age groups 

N 

1 to 10 

Number of age groups assumed to 
occur in the sample. 

Number of length groups 

N 

2 to 100 

Number of length groups in the 
simulated samples. 

L* 

N 

1 to 10 5 

Asymptotic length. 

CV ofU 

N 

0 to 99 

A measure of the assumed 
variability of L,, i.e. CV = 
o(L,)'100/L x . 

K 

N 

0.1 to 20 

Curvature parameter of the VBGF. 

CV of K 

N 

0 to 99 

A measure of the assumed 
variability of K, i.e. o(K)100/K. 

C 

N 

Oto 1 

The amplitude of the seasonal 
oscillations of the growth curve. An 
input of zero will automatically set 
the Winter Point (WP) to zero. 
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Table 5.1 (cont). Required single field inputs for simulating length-frequency data. 


Field name 

Type 

Limits 

Remarks 

CV of C 

N 

0 to 99 

A measure of the assumed 
variability of C, i.e. o(C)-100/C 

WP 

N 

Oto 1 

Time of the year when growth is 
slowest. 

CV of WP 

N 

0 to 99 

A measure of the assumed 
variability of WP, i.e. 

o(WP)100/WP. 

L50* 

N 

Oto 10 6 

50% retention length (left hand 
selection and/or recruitment). 

L75* 

N 

Oto 10* 

75% retention length (left hand 
selection and/or recruitment). 

R50* 

N 

Oto 10* 

50% retention length (right hand 
selection or de-recruitment). 

R75* 

N 

OtolO 6 

75% retention length (right hand 
selection or de-recruitment). 

F maximum* 

N 

Oto 99 

The maximum fishing mortality 
over age group that will be applied 
to the stock. 


*In migratory stock, these fields are area-specific, i.e. values are defined for 
each area the fish migrates to. 
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Table 5.2. Other required inputs when running a simulation routine (Options 1 
and 2). 


Field name 

Type 

Limits 

Remarks 

Seed value 
for randomization 

N 

Oto 100 

Use the same seed value to 
obtain exactly the same values 
as in a previous run. 

No. of samples 

N 

0 to 50 

This refers to the number of 
sampling months to simulate. 
Note that as more simulations 
are required, computing time 
increases exponentially. 

No. of months 
between samples 

N 

1 to 11 

Number of months separating 
simulated samples. 

Run time identifier 

C 

9 char. 

This will be recorded, when 
stored to a FiSAT file in 
OTHER FILE IDENTIFIERS. 
Provide inputs to this field that 
will clearly identify the 
simulation run. 

Recruitment strength 

N 

Oto 100 

The inputs for this table are in 
relative terms, i.e. the strength 
of recruitment in a given month 
vis-a-vis other months. FiSAT 
will proportionally adjust all 
values such that the largest 
entry will be equal to 1 . 

Relative strength 

N 

Oto 100 

Relative strength of each age 
group; as with the previous 
table entry, FiSAT will 
proportionally adjust all values 
such that the largest entry will 
be equal to 1 . 

Natural mortality 

N 

0.01 to 20 

The natural mortality affecting 
each age group. 
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Table 5.3. Additional inputs for running a simulation routine incorporating 
migrations. 


Field name 

Type 

Limits 

Remarks 

L50 when 
migration starts 

N 

Oto 10 6 

50% retention length (left-hand 
selection) at the start of 
migration out of the nursery 
area (Area I). 

L75 when migration starts N 

Oto 10 6 

75% retention length (left-hand 
selection) at the start of 
migration out of the nursery 
area (Area 1 ). 

Where to go from Area 1 

N 

2 to 5 

Destination (Area 1+i) at the 
start of migration out of Area 1. 

Stay in Area 1+i 

N 

1 tolO 5 

Duration of stay in Area 1+i 
(in months). 


Following the above inputs, FiSAT will prompt the user to identify the 
subsequent areas of migration (2 to 5) and durations in the specified area (1 to 
10 5 ). 


Table 5.4. Summary and sequence of procedures which defines a Monte Carlo 
simulation to generate length-frequency samples. 


Probability 

Input 

Probability' of 

Output 

distribution 

Parameters 

occurrence 


Recruitment 

no. of age groups 
(t n ) and year class 
strength (Cyj); if 
year class strength is 
not specified, class 
strength is assumed 
equal for all year 
classes 

Cyj/ZCy, 

the age 
grps. where 
the fish 
should be 
placed 

Month of 

relative monthly 

MRj/EMRj 

day is 

birth 

recruitment (MR;) 


assumed to 
be the 1 5th 
of the 
month and 
birthdate 
(tz) is 
defined 
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Table 5.4 (Cont.). Summary and sequence of procedures which defines a Monte 
Carlo simulation to generate length-frequency samples. 


Probability 

distribution 

Input 

Parameters 

Probability of 
occurrence 

Output 

Growth para- 
meters 

growth parameters 
and their corres- 
ponding standard 
deviations o(x) 

l/(x+2o(x)) 

defines the 
growth 
curve of 
a fish 

Migration 

the ages at which the 
fish is expected to 
stay in the fishing 
ground 

(n.a.) 

checks if 
the fish 
has left or 
arrived at 
the fishing 
ground 

Mortality 

Natural mortality (M t ) 
and selection data 
SEL t at age t, and the 
maximum fishing morta- 
lity the sample is 
exposed to during its 
life < F max) 

e‘ z , where 

Z-£((F max 'P() 

+M t ) 

checks if 
the fish 
survived 
to that age 
or died due 
to fishing 
or natural 
mortality 

Gear 

selection 

Probability of capture 
(i.c. selection and 
recruitment curve) by 
length or age groups 
_!Ei> 

Pt 

checks if 
the fish 
was caught 
by the gear 
or not 


Outputs Both options of this routine can 

generate time series of length 
frequencies that are saved 
automatically. 


User interface The user interface of this routine 
contains a series of tabs (Fig. 5.1) to 
encode the parameters. Default values 
are provided by FiSAT II that can be 
modified. The resulting length 
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frequencies are automatically saved for 
use in various routines of FiSAT II. 

Remarks The time required to run this simulation 

routine depends on the inputs provided, 
and may range from a few minutes to 
several days. However, it is possible to 
interrupt any simulation by pressing any 
key (warning: intermediate results are 
not saved, and the simulation must be 
restarted from scratch if interrupted). 



Fig. 5.1. User Interface when executing the simulation routine to generate length 
frequencies. 


Sample weight estimation 

This routine allows users to estimate sample weights as 
required to estimate raising factors, or for other purposes. 
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Required file 


Input parameters 


Functions 


Outputs 


User interface 


Length frequencies with constant class 
size. 

a and b coefficients of the length-weight 
relationship. 

The mean weight of the sample (Ws) is 
computed from 

w.-2* i -M l / 2 N , 

1 I 


where 


Nj is the frequency count, and w> is the 
mean weight of the fish in class i 
computed from 



where a and b are the coefficients of the 
length-weight relationship and Lj and 
L,^| are the lower limit and upper limit 
of length class i, respectively; 

A table of computed weights for each 
length class and the weight of the 
sample. 

The user interface contains two tabs 
(Fig. 5.2 and Fig. 5.3). The “General” 
tab (Fig. 5.2) identifies the file and 
provides the fields to enter the ‘a’ and 
‘b’ coefficients of the length-weight 
relationship. 
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' Sample Weight Estimation— 


Sample Weight Estimation 


*1 


General | Computed Wetfrtt) 


Data 


Fleneme 

|C SPiogiam Flei\F6AT ll\Data\Samplelq 

Species name 

|Upeneus mokjccema 


Othet idertfos 

|RAGAY GULF 1980 


Clan size: 

|1 Largest class see: |185 

Ural of measurement |cm 


Length-Wetfrl Relebonshp (W - a ■ L*b) 


coal a | 0 001 coef be | 3 


Cgmpute girt 


Close 


Fig. 5.2. User interface to encode inputs for sample weight estimation routine. 


The second tab (Fig. 5.3) is the results 
of the analysis. Note the presence of the 
“Compute” command button. After file 
identification and encoding the required 
coefficients, clicking this button will 
commence the calculations. 

Also, note the presence of the “Update 
Samples” command button. Clicking 
this command button will update the 
sub-header of the length frequency file, 
replacing whatever value is stored. This 
function should, therefore, be used with 
caution. 
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Fig. 5.3. Outputs of the sample weight estimation routine. Note that the sample 
sub-header can be automatically updated using the computed sample weights. 


Remarks The unit for the computed weights is 

determined by the coefficient a, 
provided as input to the routine. 

[The above equation provides unbiased 
estimates of w; estimating w, as the 
weight corresponds to L , leads to bias]. 

Reading Beyer (1987) 


VBGF and L/F plot 

This routine provides the option to view length- or weight- 
frequency data as a 2-dimensional plot (Fig 5.4). 

Required file Length or weight frequencies with 

constant class size. 
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Input parameter 

Function 

Outputs 

User interface 


None 

None 

A plot of the histogram in two 
dimensions. 

The user interface for this routine 
contains two tabs (as in the previous 
routine) (Fig. 5.4). The first tab 
identifies the length frequency file to 
plot and text boxes to encode the 
growth parameters to use in plotting the 
VBGF curve (if applicable). 

There are several options that the user 
may invoke in plotting the length 
frequency histograms and the growth 
curve. An option is provided to plot the 
restructured data (see Chapter 4 on data 
restructuring as employed in 
ELEFAN I). However, the option to 
plot using “Solid Bar Graphs” will be 
disabled if restructured data are to be 
plotted. 
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VBGF and Length Frequency Plot 


QoU and Graph Options j graph (Output)] 
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Species name 

|Upeneus moiuccensts 



Other (Jenifers 

fRAGAYGULF 1960 



Class tee 

1 
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Asymptotic length (loo) 
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OX 
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.J 

•Voter Point* 1 |WP) 

OX 
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Startsng length fSL| 

90 


Plot VBGF | frw | Clow 


Fig. 5.4. User interface for encoding parameters for the VBGF plot and selecting plotting 
options. 


The resulting graph (Fig. 5.5) may be 
printed or saved to file (from the Print 
command). The checkbox in the lower 
left comer of the graph, “Display all” is 
enabled only if more than one year of 
data has been detected. In this case, 
users may use the scroll bar to view 
other parts of the resulting plot or view 
the complete file in one screen by 
checking the “Display all” checkbox. 
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Fig. 5.5. Resulting VBGF plot overlaid by the length frequency histograms. 


Remarks Users should rely on this routine to 

assess visually the progression of modes 
in a series of samples before attempting 
to estimate growth parameters from 
such series. 


Maximum length estimation 

This routine estimates the (un-sampled) maximum length 
of fish (L max ) in a population, based on the assumption that 
the observed maximum length of a time series of samples 
does not refer to a fixed quantity but rather, represents a 
random variable which follows a probabilistic law. 

Required file Time series of length frequencies with 

constant class size (the routine only 
uses the extreme length of each 
sample). 

Input parameter None 
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Functions 


Output 


User interface 


L’max > s estimated from a set of n 
extreme values (L*, the largest 
specimen in each sample of a file) using 
the (Type I) regression: 

L* = a + 1/a ■ P 

where P is the probability associated 
with the occurrence of an extreme 
value, 1/a a measure of dispersion, and 
L’max * s the intercept of the regression 
line with the probability associated with 
the n th observation (note that the scale 
used for P is non-linear, i.e. corresponds 
to that used for extreme value 
probability paper). 

P is computed for any extreme value 
from P = m/(n+l) (Gumbel 1954) 
where m is the position of the value, 
ranked in ascending order and n is the 
number of L* values. 

Predicted L’ max for a given sample size 
n, and its confidence limits (here for 
95% probability). 

The user interface of this routine 
contains three tabs (Fig. 5.6 and 
Fig. 5.7). The first tab identifies the file 
to use in the analysis and summary of 
the results. The second tab displays the 
probability plot with the corresponding 
confidence region. The third tab 
summarizes the result in table format. 
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Remarks This routine may also be invoked to 

estimate L’ max for inputs to routines 
described earlier (e.g. Ault and 

Ehrhardt's model) or as an initial 
estimate of L,. 

Readings Formacion et al. (1991) 



Fig. 5.6. User interface used for the prediction of the maximum length from 
extreme values. This tab identifies the file and summarizes the results. 
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^Maximum Length Estimation... 


Piediction of the Maximum Length from Extieme Values 


General Plot j Values] 


Extreme lengths (cm) 

D — * KJ CO -U 

O O O O 














































001 .1 .5 .9 96 99 999 

Cumulative probability 


Compute | Print.. | Close 


Fig. 5.7. Facsimile representation of the resulting analysis of extreme values. 


Growth performance indices 

This utility provides the facility to rapidly compute the 
growth performance index (1) a (if asymptotic weight is 
given) or (2) o’ (if asymptotic length is given). 

Required file None 

Input parameters Option 1 : o 
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Functions 


Outputs 
User interface 


Growth parameters W, (asymptotic 
weight; g) and K (VBGF curvature 
parameter; year 1 ) 

Option 2: o' 

Growth parameters L« (asymptotic 
length in cm) and K (VBGF curvature 
parameter; year 1 ) 

Option 1: o 

0 = log lo (K) + 2/3-log lo (W,) 

Option 2: o' 

o' = log 10 (K) + 21og 10 (L S[ ) 
o and/or o' values. 

The user interface of this routine 
(Fig. 5.8) provides fields for the user to 
encode known parameters. Depending 
on what is available, the other 
parameters will be computed. For 
example, if L, and the growth constant 
K are given, phi prime (o') will be 
computed. If L, and phi prime are 
given, the growth constant K will be 
computed. The same is true with 
regards to the relation of phi, K and W„ 
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j > Growth Performance Index... 




T he parameters 0 (-log K + 2/3 log Woo) and 0 ' (-log K + 2 log Loo) have been 
shown to be remarkably constant between different populations of the same 
species, as long as similar units and definitions ate used. eg. cm and total length 
for Loo. g and live weight for Woo. 1 /year for K and base 1 0 for the logarithms 


The routine below can be used to estimate 0 or 0 '. given growth parameters, or 
K. given 0 ' and Loo. or 0 and Woo 


Growth Parameters 
Asymptotic length (Loo): 
Asymptotic weight (Woo): 
Curvature parameter (K): 


Re-compute 


Growth Performance Indices 
0 ': | 2644 

0 I T286 


Close 


21 

85 

T 


Fig. 5.8. User interface in the calculation of growth indices. 


Remarks None 

Readings Munro and Pauly ( 1 983) 

Pauly (1979) 

Pauly (1981) 

Pauly and Munro (1984) 


Regression analysis 

FiSAT II provides a routine to perform linear regression 
analyses as commonly used in fisheries stock assessment. 
Several options are provided; 

( 1 ) Y = a + b-X (Type I, or AM regression) 

(2) Y = a' + b’-X (Type II, GM regression) 

(3) log(Y) = a + b-log(X) 

(4) ln(Y) = a + b-X 

(5) Y = a + b-ln(X) 

(6) Y/X = a + b-X 

(7) ln(( 1 -Y)/C) = a + b-X 
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Required file Two-column regression file with at least 

three pairs of data. 

Input parameter Options 1 to 6 : None 

Options 7 to 9 : A constant C (>0 to 
999). 

Functions As defined by the option labels, where a 

is the Y-intercept, b is the slope, and C 
is a constant. 


Outputs 


An XY-plot of the points, and the 
regression line (Fig. 5.9). 


xj 




£rrt 


£kne 


Fig. 5.9. User interface in regression analysis of a XY file. 


User interface This routine requires an active XY 
Regression file (*.xyf). The options 
presented above can be accessed from 
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Remarks 


the dropdown list. The results (standard 
regression statistics) are given on the 
right hand side of the interface. 

The options available in this routine 
may be applied to various scenarios. 
Option 3, for example, may be used to 
establish the length-weight relationship 
of fish; Option 4 can be used for catch 
curve analysis, where ln(N) is plotted 
against absolute or relative age; 
Option 6 is typically used for surplus 
production models whose parameters 
are often estimated through a plot of 
catch per unit effort (C/f) against effort 
(f); Option 7 may be used to linearize 
complex functions, such as the VBGF 
in which case, Y = L t , C = L, and X = t. 
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APPENDIX A. Terms and variables 


AL 


AL/At 

At 


A6 

6 

6 2 

$ 

2 

ii 

0 

0 * 


a 

a' 

AM 

ASP 


B 


length increment; or 

width of length class in grouped data; or 

difference between two successive mean 

lengths 

growth rate expressed as difference 

time difference, e.g., the time needed by an 

average fish to grow from the lower to the 

upper limit of a length class 

difference between two successive standard 

deviations 

standard deviation of variates, used as a 

measure of their dispersion 

variance; the square of 6 

US dollars, or any other monetary unit 

summation sign 

pi = 3.1415... 

phi, i.e. a weight-based index of growth 
performance (0 = log ]0 (K) + 2/3 logi 0 (W,)) 
phi-prime, i.e. a length-based index of 
growth performance ( 0 ’ = log 10 (K) + 
21og 10 (L x )) 

Y-intercept in a Type I, or AM linear 
regression, or 

multiplicative term in a length/weight 
relationship 

Y-intercept in a Type II, or GM linear 
regression 

arithmetic mean; used to characterize "Type 
I" regressions 

"Available Sum of Peaks"; the sum of 
available "points" in a file restructured for 
analysis with the ELEFAN I routine 

biomass, or stock size in weight 
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b 

b' 

B/R 

B'/R 


exponent of a length-weight relationship, or 
slope of Type I (AM) linear regression 
slope of a Type II (GM) linear regression 
biomass per recruit 
relative biomass per recruit 


C 


c.v. 

C/f 
cl, c2 

CGA 

Ch 

C i,A 

C i,B 

cm 

C, 


refers to fields that will accept alpha- 
numeric entries when used to describe the 
characteristics of a field, or 
catch in numbers, or 

parameter expressing the amplitude of 
seasonal growth oscillation in the VBGF, or 
a constant 
the fraction LJL* 

coefficient of variation, i.e. C.V. = 6/X (also 
expressed in %, i.e. C.V. = 6 • 100/X) 
catch per unit of effort (also: CPUE) 
multipliers for estimating Z and its standard 
error using one of Hoenig's methods 
Colour Graphic Adapter 
refers to limited choice fields (i.e. users have 
to choose from a list). 

cumulative catch in numbers from length i 
to L„ 

cumulative catch in numbers for mesh size 
m A 

cumulative catch in numbers for mesh size 
m B 

centimetre 

terminal catch, as used in VPA 


D 

DD 

DDD 


Fraunhofer diffraction function in 

Shepherd's method, or 

dimension, as in "2D", "3D" 

two digit number denoting the day of a 

month, or 

degrees latitude 

three digit number denoting the degrees 
longitude 
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d.f. 


e 


E 

E„.l 


EGA 

EPSON 

ESP 

EXP 

F 

f 

f-factor 


f! 

Ft 


degrees of freedom, i.e. the "real" number of 
cases available for testing a statistical 
hypothesis 

base of the natural (or Napierian) 

logarithms; e = 2.71828... 

exploitation rate; E = F/Z 

level of exploitation at which the marginal 

increase in yield per recruit reaches 1/10 of 

the marginal increase computed at a very 

low value of E 

exploitation level which will result in a 
reduction of the unexploited biomass by 
50% 

exploitation level which maximizes Y/R or 
Y’/R 

Enhanced Graphic Adapter 

registered trademark of Seiko Epson Corp., 

Japan 

"Explained Sum of Peaks"; the points 
"explained" (or hit) by a growth curve traced 
by the ELEFAN 1 routine 
exponent 

instantaneous rate of fishing mortality 
fishing effort 

factor used as a multiplier to simulate a 
change in effort level, for a defined fishing 
regime 

index for fleet 

terminal fishing mortality, as used in VPA 
and cohort analysis 


g - gram 

GM - geometric mean, used in Type II regression 

HD disk - High-density disk with a capacity of 

1.44MBytes for V/i" disk and of 1.2MBytes 
for 514" floppy disk. 
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HERCULES - Hercules graphic adapter for monochrome 

screens, with a resolution of 729 by 348 
pixels. 

HP - registered trademark of Hewlett-Packard 

Co., USA 

i - symbol or subscript used for counting items 

(samples, means, etc.) 

IBM - Registered trade mark of International 

Business Machines, Corp., USA 

K. - curvature parameter of the VBGF 

k - the number of parameters estimated by a 

given procedure 

L - "length" of a fish, shrimp, etc. (length itself 

is defined differently, depending on what is 
measured, see TL, SL, FL, etc.) 

L - mean length of fish, computed from L' 

upward, or 

- mean of two or more lengths, e.g., mean of 
length at tagging and at recapture 

L' - a length not smaller than the smallest length 

of fish fully represented in catch samples; 
used to compute L 

L* - largest observed specimen in a sample 

L/F - length-frequencies or length-frequency 

sample 

L 2 5 - length at which 25% of the fish will be 

vulnerable to the gear (left-hand selection) 
L 50 - length at which 50% of the fish will be 

vulnerable to the gear (left-hand selection) 
L 75 - length at which 75% of the fish will be 

vulnerable to the gear (left-hand selection) 
LAN - Local Area Network 

L c - mean length of fish at first capture; 

equivalent to L 50 

L m - length at tagging or marking 

L mass - mean length at first maturity (or "massive 

maturation") 
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L 


max 


In 

log 

L r 

V 

La 

l b 

L t 

L. 


M 

m 


m A 

m B 

ML 

MM 


MPA 
MS DOS 


n 

N 


maximum length reached by the fish of a 

given stock; may also be predicted from the 

largest specimens of several samples using 

the extreme value theorem 

smallest length represented in one or several 

samples 

log e , logarithm of base e 

log jo, logarithm of base 10 

length at recapture, or 

mean length at first recruitment 

computed length at recapture given growth 

parameters (L, and K) and length at marking 

optimum length for mesh size m A 

optimum length for mesh size m B 

(mean) length at age t 

asymptotic length, i.e. the (mean) length the 
fish of a given stock would reach if they 
were to grow forever 

instantaneous rate of natural mortality, i.e. 
due to all causes except fishing 
mesh size, or 
metre 

gillnet mesh size 

another gillnet mesh size, with m B >m A 
"mid-length" or length class midpoint 
two digit number denoting the month of a 
year, or 

minutes in latitude and longitude 
modal class progression analysis 
disk operating system for IBM PCs or its 
compatibles, or 

registered trademark of Microsoft Corp., 
USA 

number of items in a sample, number of 
cases investigated, etc. 
number of fish in a given size class of a 
catch sample, or 
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N, 


refers to numeric fields, i.e. indicates that 
only numbers can be entered; 
number of fish in the oldest age group of a 
cohort or population ("terminal population") 


P 

Pl 

Pi 


PC 

Prompt 


probability of capture or occurrence on the 
fishing ground 

probability of capture for length (or mid- 
length) L 

first point of a length-converted catch curve 
included in the computation of Z; this point 
is by definition the first where the 
probability of capture is 1 
Personal Computer; also microcomputers 
a software message or signal inviting the 
user to enter data, or to perform an operation 


R 


50 


R 


75 


RAM 


Rn 


product-moment correlation coefficient 

coefficient of determination 

length at which 50% of the fish will no 

longer be vulnerable to the gear (right-hand 

selection, or deselection) 

length at which 75% of the fish will no 

longer be vulnerable to the gear (right-hand 

selection, or deselection) 

Random Access Memory; a part of the 
memory of the computer where the program 
and the data are loaded. 

"goodness-of-fit" index of the ELEFAN I 
routine (=10 ESp/ASP /10) 


S 

SF 

SI 

SL 


SS 


score function in Shepherd's method 
selection factor 
separation index 

starting length; one of the two coordinates 
used to locate a growth curve in the 
ELEFAN I routine 

starting sample; the other coordinate used to 
locate a growth curve in the ELEFAN I 
routine. Jointly, SL and SS define the 
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s.e. 

Sj, S2, S3, S4 

s 

J max 

SSE 

t 


T 

t' 


tm 

*-tnass 

tjnax 


tr 

's 

t z 


v i 

VBGF 

VGA 


location of a pre-selected point of a growth 
curve, or 

seconds in latitude and longitude 

standard error of a statistic 

variables used for estimating the probability 

of capture under the logistic model 

maximum score in Shepherd's method for a 

range of L, and K parameters 

sum of squared errors; a measure of 

dispersal from the mean 

a given time or age (normally expressed in 
years) , or 

absolute age of a fish, e.g., as estimated 

from daily otolith rings, or 

age corresponding to L t 

mean annual habitat temperature, in °C 

relative age of a fish, defined as t' = t - t 0 

mean age at first capture, corresponding to 

Lc 

mean age at length i 
age at marking, corresponding to L m 
mean age at massive (*• first) maturity 
longevity (in the wild) 

the "age" fish would have had at length zero 
if they had always grown according to the 
VBGF; t 0 generally has a negative value, but 
does not usually express "prenatal growth" 
mean age at recruitment 
parameter of the seasonally oscillating 
version of the VBGF (see WP) 
in Shepherd's method: origin of the VBGF 
in calendar time, (expressed as fraction of a 
year); here replaced (without affecting other 
results) by a starting point, defined by SS 
and SL 

total of estimated value at length i 
von Bertalanffy Growth Function, either in 
original or seasonally oscillating form 
Virtual Graphic Array 
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VPA 

- Virtual Population Analysis 

W 

- mean weight of fish in catch samples, 
computed from W' upward 

w 

- mean weight of fish within a given length 
class 

WF 

- weighting factor assigned to an observation 

WP 

- "Winter Point"; in the scasonalized VBGF, 
the time of the year when growth rate is 
slowest; equivalent to t s + 0.5 year 

WC 

- total weight of the catch 

WS 

- total weight of the sample 

W, 

- asymptotic weight, i.e. the (mean) weight 
the fish of a given stock would reach if they 
were to grow forever 

X 

- any variable (often used for the abscissa in 
2-dimensional plots) 

X 2 

- chi-square statistics 

Y 

- yield, catch in weight 

Yj 

- yield at length i 

y 

- any variable (often used for the ordinate in 
2-dimensional plots) 

YY 

- two-digit number denoting the year 

Y/R 

- yield per recruit 

Y’/R 

- relative yield per recruit 

Y/R max 

- maximum yield per recruit achievable under 
a given fishing regime 

Z 

- instantaneous rate of total mortality 
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